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We revisit the problem of dynamical response in spin-charge separated one dimensional quantum 
fluids. In the framework of Luttinger liquid theory, the dynamical response is formulated in terms of 
noninteracting bosonic collective excitations carrying either charge or spin. We argue that, as a result 
of spectral nonlinearity, long-lived excitations are best understood in terms of generally strongly 
interacting fermionic holons and spinons. This has far reaching ramifications for the construction 
of mobile impurity models used to determine threshold singularities in response functions. We 
formulate and solve the appropriate mobile impurity model describing the spinon threshold in the 
single-particle Green’s function. Our formulation further raises the question whether it is possible 
to realize a model of noninteracting fermionic holons and spinons in microscopic lattice models of 
interacting spinful fermions. We investigate this issue in some detail by means of density matrix 
renormalization group (DMRG) computations. 
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I. INTRODUCTION 

Understanding the essential features of a quantum 
many-body system usually entails finding a simple expla¬ 
nation of its low-energy spectrum in terms of weakly in¬ 
teracting, long-lived quasiparticles. For example, in Lan¬ 
dau’s Fermi liquid theory the quasiparticles are fermions 
carrying the same quantum numbers as an electron, 
namely charge e and spin 1/2. However, it is well estab¬ 
lished that the long lived excitations in strongly corre¬ 
lated systems may carry only a fraction of the quantum 
numbers of the elementary constituents. In fact, frac- 
tionalization is often invoked as a route towards exotic 
phases of matter such as spin liquids or high-temperature 
superconductors^. 

Perhaps the most prominent example of fractionaliza- 
tion is spin-charge separation in one dimensional (ID) 
quantum fluids known as Luttinger liquids^. The hall¬ 
mark of these theories is a low energy spectrum described 
by two decoupled free bosonic fields associated with col¬ 
lective spin and charge degrees of freedom, respectively. 
On the experimental side, the most direct evidence for 
spin-charge separation involves the observation of mul¬ 
tiple peaks associated with spin and charge collective 
modes in dynamical response functions at fairly high 
energies^^®, beyond the regime where Luttinger liquid 
theory is applicable. Spin-charge separation is known 
to persist at high energies for integrable theories such as 
the ID Hubbard^^^® and 1/r^ models. In these 

cases any eigenstate with a finite energy in the thermo¬ 
dynamic limit can be classified in terms of elementary 
excitations called holons (which carry charge e and spin 
0) and spinons (which carry charge 0 and spin 1/2), along 
with their bound states. A convenient way for describing 
such excitations as well as their scattering is to define 
corresponding creation and annihilation operators Z\{9) 


and Za{0), where a labels the different types of elemen¬ 
tary excitations and 0 is a variable that parameterizes 
the momentum Pa{0) (the form of this function depends 
on the particular model under consideration). As a con¬ 
sequence of integrability, creation and annihilation oper¬ 
ators fulfil the Faddeev-Zamolodchikov algebra^®’^® 

Za{e^)Zb{e2) = sii{ei,e2)Zd{e2)z,{ei), 
Za{9i)zl{B2) = 2n5{ei - e2)5a,b 

+st{e2,e,)z\{e2)z,{e,). (i) 

Here 5'/)((0i,02) is the purely elastic two-particle S- 
matrix. The corresponding elementary excitations are 
infinitely long lived, but generally strongly interacting, 
as can be seen from their scattering matrices®^^®. More¬ 
over, their quantum numbers, e.g. charge 0 and spin 1/2 
for spinons in the Hubbard model, differ from those of 
the collective bosonic spin modes in the Luttinger liquid 
description. It is natural to assume that breaking inte¬ 
grability slightly will not generically lead to a qualitative 
change in the nature of elementary holon and spinon ex¬ 
citations, but merely render the lifetimes finite. 

How to reconcile this picture emerging from the ex¬ 
act solution with Luttinger liquid theory? In the latter, 
the nature of elementary excitations is obscured by the 
fact that, due to the linear dispersion approximation, the 
spectrum is highly degenerate and allows for many inter¬ 
pretations, which ultimately all give the same results for 
physical observables. Examples are chiral Luttinger liq¬ 
uid descriptions of quantum Hall edges^^ and the low- 
energy excitations of the Hubbard model. The latter 
can be understood both in terms of interacting, fermionic 
holons and spinons^’^^^^^, and in terms of noninteracting 
bosons associated with collective spin and charge degrees 
of freedom^®^^^. 

Going beyond the linear dispersion approximation is 
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expected to remove ambiguities in the quasiparticle de¬ 
scription of the spectrum: there will be a particular 
choice that maximizes the lifetimes of elementary exci¬ 
tations. In integrable models such as the Hubbard chain 
these lifetimes are infinite. 

Over the last decade it has been established in a series 
of works that going beyond the linear dispersion approx¬ 
imation is essential for correctly describing the dynamics 
of one dimensional models with gapless excitations^®^®®. 
This has resulted in the so-called “nonlinear Luttinger liq¬ 
uid” (nLL) approach to gapless ID quantum liquids; see 
Refs. 46 and 54 for recent reviews. The basic reason for 
the failure of linear Luttinger liquid theory is that at fi¬ 
nite energies the running coupling constants of irrelevant 
perturbations such as band curvature terms are in fact 
different from zero. Taking them into account in pertur¬ 
bation theory leads to infrared singularities. These need 
to be resummed to all orders in the coupling constants, 
which gives rise to new, momentum-dependent exponents 
in response functions. 

A key ingredient of the nLL method is the identifi¬ 
cation of quasiparticles describing excited states both at 
high and at low energies. This is straightforward for spin¬ 
less fermions, because the quasiparticles are adiabatically 
connected to free fermions^®. The spinful case is consid¬ 
erably more involved due to the onset of spin-charge sep¬ 
aration for arbitrarily weak interactions and the concomi¬ 
tant qualitative change in the nature of the elementary 
excitations compared to the noninteracting limit. It was 
realized in Refs. 49-52 that in many important cases ex¬ 
act results can be obtained by using a phenomenological 
model of weakly interacting fermionic holons and spinons, 
whose dispersions delineate the edges of the support of 
the dynamical response function under consideration. By 
construction these excitations are different from the true 
elementary excitations in the Hubbard model. In partic¬ 
ular, they carry different spin and charge quantum num¬ 
bers. It is then an obvious but crucial question how to 
reconcile this approach with the exact solution of inte¬ 
grable models like the Hubbard chain. 

In this work we develop a new approach to deriving 
mobile impurity models for studying dynamical corre¬ 
lations in gapless models of spinful fermions. We first 
carry out a direct construction of the “physical” holon 
and spinon fields in the limit of weak electron-electron 
interactions. This results in a representation of spinful 
nonlinear Luttinger liquids in terms of strongly interact¬ 
ing holons and spinons, which is in direct accord with 
known results obtained in integrable cases. Recalling 
that the spinless fermion case was best understood by 
considering weak interactions, we address the construc¬ 
tion of a microscopic model of interacting electrons giv¬ 
ing rise to a theory of noninteracting fermionic holons 
and spinons at low energies, but in presence of spectral 
nonlinearities. We then derive mobile impurity models 
to analyze dynamical response functions in our formu¬ 
lation. We demonstrate explicitly how to recover results 
obtained previously by means of the approach of Refs. 51 


and 52. Finally, we address the question how to realize 
a lattice model of interacting electrons that gives rise to 
noninteracting holons and spinons at low energies. 

II. SPINLESS FERMIONS 

Before we approach the problem of defining quasipar¬ 
ticles for spin-1/2 fermions, it is instructive to review the 
case of spinless fermions. For concreteness, consider the 
simplest lattice model of interacting spinless fermions: 

L 

H = -f h.c.) + V'^njUj+i. (2) 

i=i j 

Here Cj is the annihilation operator for a fermion at 
site j, Uj = cJjCj is the number operator, t is the hop¬ 
ping parameter, and V is the nearest-neighbor interaction 
strength. This model has a U(l) symmetry, Cj —)• 
a G K, associated with conservation of the total number 
of fermions N. Moreover, the model is integrable and 
the exact spectrum for arbitrary values of V can be cal¬ 
culated from the Bethe ansatz solution®^’®®. There is a 
gapless phase extending between —2<V <2. AtV = 2 
the model is equivalent to the SU(2)-symmetric Heisen¬ 
berg spin chain®^. 

A useful starting point for analytical approximations 
is the noninteracting model with V = 0. In this case, the 
elementary excitations are free fermions with dispersion 
relation €o(k) = —2tcosk — ^ measured from the Fermi 
energy /i. Low-energy excitations are particles and holes 
with momentum close to the Fermi points, k « ±kp. 
The Fermi momentum is related to the average density 
by kp = TrN/{Lao), where uq is the lattice spacing. 

At weak coupling, Vjt 1, standard bosonization can 
be used to derive an effective low-energy theory for model 
(2), see e.g. Refs. 25, 26, and 59. One starts by taking 
the continuum limit and projecting the fermionic field 
onto states with momentum near the Fermi points. This 
leads to the right- and left-moving components in the 
mode expansion 

Cj -k L{x)]. (3) 

The effective Hamiltonian in terms of R and L fermions 
reads 

R = J dx[—v'R^idxR + v'L^idxL 

+gR^RL'^L + R;„{x)], (4) 

where all operators are normal ordered with respect to 
the noninteracting Dirac sea. To first order in V, we 
have the parameters v' ~ 2tao sinkp (l -I- sinfe^;’) and 
g « 4:Vao sin^ kp. The coupling constant g is the only 
marginal interaction. The term Ri„{x) contains nonlin¬ 
ear dispersion terms and other interactions which are ir¬ 
relevant in the renormalization group (RG) sense. In the 
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standard Luttinger liquid approach this term is dropped, 
which corresponds to linearizing the dispersion around 
the Fermi points 

The bosonization formula for chiral spinless fermions 
reads 

= (5) 

= ( 6 ) 

where 77 , fj are Majorana fermions and ip{x), (p{x) are chi¬ 
ral bosons that obey the commutation relations 

[(p{x),<f{y)] = 0, (7) 

[ip{x),Lp{y)\ = 27risgn(a; - y) = - [ip{x), i^{y)]. ( 8 ) 

Throughout this paper we use “CFT normalizations” for 
bosonic vertex operators 

^-5--^. (9) 

[x - yY°‘ 

A consequence of employing these conventions is that 
vertex operators are dimensionful 


Thus, an elementary excitation with charge g = 1 corre¬ 
sponds to a kink of amplitude 7 r -\/8 in the bosonic field 
$(x). As is well known, the model in Eq. (15) correctly 
captures the long-distance asymptotic decay of correla¬ 
tion functions for any ID system in the Luttinger liquid 
universality class^^^^^. 

The limitations of Luttinger liquid theory appear when 
one considers dynamical response functions at small but 
finite frequency and momentum. At finite energy scales 
it becomes necessary to take the irrelevant perturbations 
"Hirr to the Luttinger model into account. At weak cou¬ 
pling and in the absence of particle-hole symmetry (i.e. 
away from half-filling in the lattice model), the leading 
corrections in (4) are the dimension-three operators 

H,,. = -^{R^dlR + L^dlL) 

RLhd^L - L^^LRYO^R + h.c.). (18) 

Here we have introduced rh~^ « m~^ H—^ sin 2kF, with 
m~^ = 2 fagCOs/cF the inverse free fermion mass, and 
gs « Hog sin2/ci^’. Bosonizing these irrelevant terms we 
obtain cubic boson-boson interactions 


dim = length""'. ( 10 ) 

We also define the canonical bosonic field 4) (a;) and its 
dual 0 (a;) by 

^{x) = ip{x) + If{x), (11) 

e{x) = ip{x) - Ifix), (12) 

which obey 

[$(a;), 0(a:O] = 47 r 7 Sgn(a: — x'). (13) 

Bosonization of Eq. (4) then leads to the Hamiltonian 

R = J dx['HLL{x)+'Hi„{x)]- (14) 

The first term. 


R-hhix) 


V 

IGtt 


Kid^eY 




(15) 


is the Luttinger model. The velocity v and the Luttinger 
parameter K are given to first order in E by n « n' 
and AT « 1 — ^ sinkp. The bosonic fields describe the 
collective low-energy density mode of the quantum fluid. 
The uniform part of the density operator is related to 
4?(x) by the bosonization relation 

Q{x) = rYx)R{x) + lYx)L{x) ~- ^^dx^{x). (16) 

TTVS 

The total charge operator is given by the integral 


with A;^ 


Rirr = A+a,$[(9,$)2 + (a,0)2] 

+A 3 5,,$[(9,$)2 - (9,0)' 

T— -I- 23.1 
L m TT J o 


(19) 

_ _ _ __r j_ _ 293 ] 

While the bosonic representation allows one to take the 
marginal interaction g into account exactly, perturbation 
theory in the nonlinear boson interactions (19) suffers 
from infrared divergences^^’"^®. The latter are associated 
with the huge degeneracy of states in the linear dispersion 
approximation: all states with an arbitrary number of 
bosons moving in the same direction carrying the same 
total momentum are degenerate^^’^®. 

A way to circumvent these difficulties in analyzing the 
nonlinear bosonic theory is suggested by reverting to the 
fermionic representation (18). At the free fermion point, 
the irrelevant interaction vanishes, 53 = 0 , and one is left 
with the quadratic dispersion term with effective mass m. 
Taking the nonlinear dispersion into account in the free 
fermion model removes the degeneracy of particle-hole 
pairs that carry the same total momentum. One can then 
approach the problem from free fermions with nonlin¬ 
ear dispersion and include interactions perturbatively'®. 
This approach reveals that the most pronounced effect of 
the interactions is to give rise to power-law singularities 
at the edges of the excitation spectrum. While a complete 
analytical solution of model (4) taking into account both 
band curvature and interaction is highly nontrivial, the 
edge singularities can be described by an effective impu¬ 
rity model in analogy with the x-ray edge problem^®’®*’. 

Consider, for instance, the single-fermion spectral 
function 


q = 


p -\-00 


dx Q{x). 


(17) 


A{oj, k) = -Im Gret(w, fc), 

TT 


( 20 ) 
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where 


G^ret 


poo 

-ij 

x(V'o|{CjH_iW>c]}|V'o) 


( 21 ) 


is the retarded Green’s function, with |^o) the exact 
ground state. We can separate the negative- and positive- 
frequency parts of the spectral function: 


A{uj,k) = A^{uj,k) + Ay{u>,k). ( 22 ) 


The Lehmann representation reads: 


A^{uj,k) = 27r^ |(■^/'n|Cfc|'0o)P<5(w + - Eq), (23) 

n 

Ay{u;,k) = 27r^ |(•^/'nl4|'0o)P<5(w - E^ + Eq), (24) 

n 

where Ck annihilates a fermion with momentum k and 
\tjjn) denotes an exact eigenstate of the Hamiltonian with 
energy En- Let us focus on the negative-frequency part 
for k < kp- For free fermions, we have A^^\uj,k) = 
d{u} — eo{k)), where eo(fc) < 0 is the energy of the particle 
annihilated below the Fermi surface. When weak interac¬ 
tions are turned on, the renormalized fermion dispersion 
e(fc) becomes a threshold of the support of A^ (oj, k), such 
that a power-law singularity develops for w < e(fc). To 
describe the edge singularity for fixed k < kp, we go 
back to the mode expansion in Eq. (3) and generalize it 
to include three patches of momentum: 


cj ^ v^[e*'=^"r(a:) + (25) 

Here the “impurity field” x^(a:) creates a hole in a state 
with momentum close to k, within a subband of width 
A ^ kp — k. The low-energy Fermi fields r{x) and l{x) 
are also defined with a cutoff of order A. In the case 
q = kp — k kp, the field l{x) can be regarded as 

the projection of L{x) into a narrower subband, while 
R{x) is split into two separate subbands corresponding 
to the low-energy mode r(x) and the “high-energy” mode 
x{x). Restricting the energy window to the vicinity of 
the threshold with a single impurity, we now have to 
calculate the propagator of x(a:): 

pOO pOO 

G'ret,<(w, A:) « —z / / dx 

Jo J — oo 

x('0o|x(O,O)x'''(La:)|V'o)- (26) 


At weak coupling, i.e. as long as the four-fermion in¬ 
teraction strength V is small, we can substitute Eq. (25) 
into Eq. (2) and bosonize the low-energy fields to de¬ 
rive an effective Hamiltonian for the single hole coupled 
to low-energy collective modes. The result is the mobile 
impurity modeF® 


^imp — 


J dx^ 




+x\s - iudE}x 
+X^x[f{(l)dx^ + /(g)i9x<p]|, 


( 27 ) 


where e = —e(fc) > 0 is the energy of the “deep hole” 
excitation, u = ^ is the velocity obtained by lineariz¬ 
ing the dispersion around the centre of the impurity 
subband, and f{q) and f{q) are momentum-dependent 
impurity-boson couplings of order V. The calculation 
of the Green’s function in Eq. (21) is made possible by 
performing a unitary transformation that decouples the 
impurity from the low-energy modes: 


U — [yifi{x)+‘ytp(x)\xAEx{E 


(28) 


The bosonic fields transform as 

ip°{x) = U(p{x)U^ = ifiix) — 2 ttjC{x), 
ip°{x) = U(p{x)U^ = if{x) + 2TTjC{x), (29) 

where 

/ OO 

dy sgn{x - y)x\y)xiy)- (30) 

-OO 

The transformed impurity field is 
d{x) = Ux{x)U^ 

— -E)Cix) ^ ^02^^ 


Note that the impurity density is invariant under the 
unitary transformation, i.e. xH^)x{^) = d^{x)d{x). 

We choose the parameters 7,7 as the solution of 


where 



„)(()■ 

= (33) 


With this choice, the Hamiltonian becomes noninteract¬ 
ing: 


Ffimp — 


V 

. IGtt - 


J dx^ 

+d^ (e — iudx)d -I- ... | 


K{dx.Q°f + 


(34) 


where ... stands for irrelevant operators which are ne¬ 
glected in the impurity model (since they only intro¬ 
duce subleading corrections to edge singularities). On 
the other hand, the expression in Eq. (26) now becomes 


^+00 poo 

Gret,<(w, fc) « —i / dte^‘^* / dx{d {0,Q)d\t,x))o 

Jo J — OO 


x{F{0,0)F^t,x))o, 


(35) 


where ()o denotes the expectation value in the noninter¬ 
acting ground state It/jo) = U\'tpo) and F{x) is the string 
operator 


F(x) = 


(36) 


The correlation function in Eq. (35) can then be cal¬ 
culated by standard methods®^. The important point is 
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that the scaling dimension of the operator F{x) changes 
continuously as a function of 7 , 7 . As a result, the effec¬ 
tive impurity model predicts a power-law singularity in 
the spectral function, 

A<(a;, k) ~ 0(e(fc) — w)|e(fc) — K (37) 

Importantly, the impurity mode x(a;) in Eq. (27) car¬ 
ries charge q = 1 because it is defined from the original 
fermion Cj at the noninteracting point. This is the parti¬ 
cle that can be identified with an elementary excitation in 
the Bethe ansatz solution for the integrable model. From 
the exact S-matrix it is known that interactions between 
these elementary excitations increases as V increases. 
Particularly at the SU(2) point, V = 2, the elementary 
excitations are rather strongly interacting. By contrast, 
the transformed impurity operator d{x) = xix)F(x) car¬ 
ries a fractional charge that depends on the interaction 
strength, since the string F{x) in general does not com¬ 
mute with the charge operator in Eq. (17). 

In the low-energy limit, an alternative approach to ob¬ 
tain the edge singularity in the spectral function was put 
forward in Ref. 44. In this approach one starts by intro¬ 
ducing fermionic quasiparticles that are asymptotically 
free at low energies. In our notation, the idea is to define 
chiral bosons (j), ^ by 

^{x) = VK{(j) + 4>), (38) 

^)- (39) 

In terms of these, the Luttinger model (15) reads 

^LL(cr) = ^P.<).)2 +(5,^)2]. (40) 

The quasiparticles R{x) and L{x) are defined by 

= (41) 

Z(a;) = (42) 

The commutator with the charge operator yields 

[g,.R'^(a;)] = VKR^ix), [q,L^x)] = '/KL\x). (43) 

Thus, the quasiparticles carry charge ^/K- On the 
other hand, this refermionization procedure removes the 
marginal interaction between the quasiparticles since the 

chiral modes are decoupled in Eq. (40)^®44 leading 

interaction is then represented by the irrelevant operator 
53 in Eq. (18), which can be neglected as a first approxi¬ 
mation in the low-energy limit. The relation between the 
original right-moving fermion and the new quasiparticle 
is 

R{x) = R{x)Fo{x), (44) 

where Fq{x) is the limit q = kp — k ^ 0 oi the string oper¬ 
ator in Eq. (36). At this point a noninteracting impurity 


mode x{x) can be introduced by projecting the free field 
R{x) into low-energy and high-energy subbands. This 
leads to a universal result for the exponent in the vicin¬ 
ity of threshold, |w — e(fc)| ^ k'^jm, which corresponds 

to Eq. (37) with parameters 7 = 

and 7 = ^ {2FK ~ result differs from the 

prediction of the linear theory 

In summary, there are two possible paths towards cal¬ 
culating edge exponents in the nLL theory for spinless 
fermions: (i) starting with free fermions, one defines low- 
energy and impurity subbands, and then turns on interac¬ 
tions between the elementary excitations in the impurity 
model; after that, the interaction with the impurity is re¬ 
moved by a unitary transformation, which introduces the 
string operator in the correlation function; or (ii) start¬ 
ing from the Luttinger model for interacting fermions, 
one refermionizes into weakly interacting quasiparticles, 
which differ from the original fermions by a string opera¬ 
tor, and then projects the quasiparticles into low-energy 
and impurity subbands. The projection onto the impu¬ 
rity model is well controlled in both paths because the 
model of interacting spinless fermions is smoothly con¬ 
nected with the free model; i.e., the parameters 7 , 7 , 
which quantify the scattering between high-energy and 
low-energy particles, vanish continuously as E —>■ 0 . 
However, it is important to emphasize the difference be¬ 
tween the original fermions, which carry unit charge of 
the U(l) symmetry, and the quasiparticles with fractional 
charge. While the latter are always weakly interacting in 
the low-energy limit, the fermions that carry the correct 
quantum numbers become strongly interacting even at 
low energies as V increases. 

Once the low-energy, weak-coupling regime is well un¬ 
derstood, the impurity model of nLL theory can be ex¬ 
tended phenomenologically to high energies, strong inter¬ 
actions and thresholds with more than one impurity, as 
has indeed been done successfully for spinless fermions®"^. 


III. SPINFUL FERMIONS 

As we have seen above, the spinless case is most easily 
understood by considering the vicinity of noninteracting 
fermions. The situation is very different in the spinful 
case. In order to understand this point in some detail, 
let us consider the particular example of the ID Hubbard 
model 

L 

iJnub = 

f = l O’ j 

(45) 

where Cj^a annihilates a fermion with spin cr =t, i at 
site j, rij^a = Cj ^Cj ^ is the number operator, and U > 
0 is the strength of the on-site repulsion. We work at 
fixed density below half-filling with zero magnetization, 
= Ni < L/2. In this case the model has a global 
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U(1)(8)SU(2) symmetry. Let us focus first on the limit 
of weak interactions U/t 1 and low energies. It is 
well known^ that the low-energy degrees of freedom are 
collective spin and charge modes respectively, i.e. 

^^Hub -'^charge + ’^spin + • ■ ■ (46) 

where the dots denote additional terms that are irrele¬ 
vant in the renormalization group sense. Crucially, as 
^^Hub is spin rotationally symmetric, "Hspin must exhibit 
a spin SU(2) symmetry. In order to parallel our analysis 
in the spinless case we wish to express Uspin in terms 
of fermionic fields carrying spin quantum numbers ± 1 / 2 . 
This is possible in an SU(2)-symmetric way only if the 
fermions are strongly interacting, i.e., the situation is sim¬ 
ilar to the V = 2 case for spinless fermions. In the charge 
sector the situation is analogous unless we work at very 
low electron densities. In order to generalize the mobile 
impurity model construction reviewed in section II to the 
spinful case, we therefore cannot work with weakly inter¬ 
acting spinful fermions, but require a model that gives 
rise to noninteracting fermions describing the collective 
spin and charge degrees of freedom. As such a model is 
not known, we proceed along the lines sketched in Fig. 1. 


weakly interacting 
spinful fermions 

I'i 


refermionize 


strongly interacting 
spin/charge fermions 


mobile impurity model 
for weakly interacting 
spinful fermions 


increase 

interactions 


t 


mobile impurity 
model 


t 


decrease 

interactions 


noninteracting 

spin/charge 

fermions 


Figure 1. Strategy map for calculating edge exponents in 
nonlinear Luttinger liquid theory for spinful fermions in the 
low-energy limit. See the main text for details. 


1. Starting with weakly interacting spinful fermions 
at low energies, we derive the corresponding model 
of strongly interacting spin and charge fermions. 

2. We then decrease the interactions in the spin and 
charge fermion model, and derive a low-energy ef¬ 
fective Hamiltonian in the vicinity of the “Luther- 
Emery point”where the spin/charge fermions be¬ 
come noninteracting. 

3. Having completed this construction, we are in a 
position to construct mobile impurity models by 
following the logic employed in the spinless case. 


4. Having constructed a suitable mobile impurity 
model, we may calculate threshold exponents by 
standard methods. 

5. Through an appropriate tuning of the parameters 
defining our mobile impurity model, we may ana¬ 
lyze the case of weakly interacting SU(2)-invariant 
spinful fermions. This is analogous to the analysis 
of strongly interacting spinless fermions with SU (2) 
symmetry based on a mobile impurity model for¬ 
mulated at weak coupling. 

A key ingredient in our approach is our “Luther-Emery 
model” of noninteracting holons and spinons. An obvious 
question is what such a theory might look like in terms 
of interacting spinful fermions. We address this issue in 
section VH. 


A. Bosonization at weak coupling 

As our point of departure we choose a general extended 
Hubbard model below half-filling, where we allow fairly 
general electron-electron interactions in addition to (45), 
provided that they are invariant under the following sym¬ 
metries: 

• U(l)(8) U(l) transformations in the charge and spin 

sectors t ot^ G ffi; 

• spin reflection 

• site parity cj^a- —t 

• translations Cj^a- Cj+i^o-. 

The kind of lattice model we have in mind is of the form 

H = iLpub ± EE 

r>l j 

+ E E ■ s,+, ± , (47) 

r>l j 

where the coupling constants {W, Jr, J/} must be such 
that the model remains in a spin-charge-separated quan¬ 
tum critical phase. Lattice models of this type can be 
bosonized by standard methods^®’^^. The generalization 
of Haldane’s bosonization formulas®^ to the spinful case 
is 

r ~ A JkFx{l-4n-2m) 

Here = 1 = — s^, oq is a short-distance cutoff, 
kp = T^N^/iLcio) is the Fermi momentum, An,m are non- 

universal amplitudes and = ’7cr(?7o-^o-»7g'^g')"'(»?rTj7s')'" 
are Klein factors (and we use notations where e.g. f ='0 
that ensure the correct anti-commutation relations. The 
bosonic fields 

^a{x) = (Pa{x) + ipa{x), (49) 

Qaix) = ipa{x) - (paix), ( 50 ) 
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with a = c, s, obey commutation relations 

0a'(x')] = 47ri(5Q,a'Sgn(x - x'). (51) 


We note that in the CFT normalizations (9) the am¬ 
plitudes An^m are dimensionful, i.e. are proportional to 
appropriate powers of the lattice spacing qq. 

In the spin-charge separated Luttinger liquid phase the 
low-energy effective Hamiltonian for extended Hubbard 
models of the type (47) is 


77 = y dx [77 ll(x) -I- 77irr(x)] , 


(52) 


Hll(x)= ^ 


77irr(x) = Ai cos + X2dx^c COS 


+ X4dx^sdxQsdxQc + • • ■ 


(54) 


Here Va are the velocities of the collective charge and spin 
modes and the corresponding Luttinger parameters. 
The contributions 77irr are irrelevant in the renormaliza¬ 
tion group sense. A complete list of irrelevant operators 
with scaling dimensions of at most four (for Kg = 1 ) is 
given in Appendix A. The velocities Va and Luttinger 
parameters Ka can be calculated exactly for the Hub¬ 
bard model, but Eq. (54) is generic for spinful Luttinger 
liquids if we regard Va and as phenomenological pa¬ 
rameters. We note that, as a consequence of spin reflec¬ 
tion symmetry, marginal interactions coupling spin and 
charge such as 


dx^s{x)dx^c{x) (55) 


are not allowed. Hence the collective degrees of free¬ 
dom at low energies are described in terms of pure spin 
and pure charge modes, rather than linear combinations 
thereof (which would be the case in presence of a mag¬ 
netic held, see e.g. Ref. 63). 


B. Refermionizing in terms of spin and charge 
fields 


The next step is to refermionize (54) in terms of spin 
and charge fermion fields. In order to see how this should 
be done, we consider the limit of vanishing interactions. 
Here the bosonization formulas simplify to (cr =t,),) 


C(7 

7?<t(x) 
Act (x) 




(56) 


where = 1. The idea is to decompose the right- 

moving spin-up electron into a right-moving holon held 
Rc and a right-moving spinon held Rg in the form 


R^{x) ~ i?c(x)e* (57) 


In a pure Luttinger liquid there are infinitely many ac¬ 
ceptable choices for R^ and string operators in Eq. (57). 
In Refs. 51 and 52, the fermions R^ are chosen so as to 
have scaling dimension 1 / 2 , in analogy with the proce¬ 
dure in the spinless case, see Eqs. (41) and (42). This 
choice is such that the particles are asymptotically free 
at low energies. However, they then carry fractional 
spin and charge quantum numbers. Such a choice is 
not the most natural one for our purposes: it is known 
that the scaling limit of the Hubbard model is given 
by the U(l) Thirring model (SU(2) at half-Ailing), see 
e.g. Ref. 26. The U(l) Thirring model is integrable, and 
the elementary excitations are known to be strongly in¬ 
teracting fermionic spinless holons and neutral spinons 
(with a known S-matrix) carrying charge =Fe and spin 
±1/2 respectively. The principle guiding our construc¬ 
tion is that charge and spin fermions created by R\ and 
R\ should carry the same quantum numbers as the el¬ 
ementary holon and spin excitations. The U(l) charges 
corresponding to these quantum numbers are 

/ + 00 

dx Qaix), (58) 

-OO 

where 

Qc{x)= [rI{x)R^{x) + lI{x)L^{x)] , 

Qs{x) = Y Sa [R^(x)i?„(x) ± L^(x)L„(x)] . (59) 

As usual these expressions are to be understood in terms 
of a standard point splitting and normal ordering pre¬ 
scription. We now require 

Qc{x) = RI{x)Rc{x) ± lI{x)Lc{x) , 

Qg{x) = Rl{x)Rg{x) + Ll{x)Lg{x), (60) 

which ensure that the spin and charge fermions have the 
desired U(l) charges 

[Qa,Riix)] = R1(x), [<?a,Ll(x)] = L1 (x). (61) 

Our refermionization prescription then reads 
7?^(x) -- r]^ Oc{x) Og{x) , 

■^ 1 ( 2 ;) ^ Vi C>c(x) Ol{x) , 

0 „(x) ~ R„(x) . (62) 

Analogous relations hold for left-moving fermions. One 
issue that arises here is that ng{x) = J^^dx'Qg{x') is 
the number of spinons on the interval [—oo, x], and there¬ 
fore string operators of the form exp [iQ!ns(x)] should be 
27r-periodic functions of a. When bosonizing string op¬ 
erators naively this periodicity is lost. A simple way of 
dealing with this issue is via the replacement®^ 

exp [iang{x)] —!> 'Y^ 6 xp [i(27rm ± a)ng{x)]. (63) 








The operators Oa{x) fulfill braiding relations for x ^ y: 
0^{x)0o.{y) = (64) 

The low-energy effective Hamiltonian (54) is expressed in 
terms of our fermionic charge and spin fields as 

% = j dx\Hc{x)+'Hs{x) A-l-Lcsix)], 

He = Rl{-iv'^d^ - ydl)R^ -f Ll{iv'^d^ - ydl)L^ 

d-gcfiR^c^c^\Lc , 

Us = Rl{-iv',d^ + iC,dl)R, + Ll{iv',d^ - iC,dl)L^ 

+gs,oR\RsL\L^ -\- gs,i{R\L^ -\- L\Rg) -f ..., 

= gi[RlRc + lIlMl, + h.c.) + .... (65) 

A crucial feature of this expression is that the coupling 
constants of the marginal interactions, 

ga,o ~ 27rwQ, -, (66) 

are 0(1) at weak coupling —t 1- Moreover gg is al¬ 

ways large as long as the spin SU(2) symmetry is unbro¬ 
ken, as in this case the Luttinger parameter is fixed at 
Kg = 1. This implies that the spin and charge fermions 
are strongly interacting. This is consistent with known 
results for the exact S-matrix of the Hubbard model*“^°. 
Moreover, the spin sector of (65) describes a massive 
Thirring model perturbed by irrelevant operators, which 
is precisely what one would expect on the basis of the 
known S-matrices for the Hubbard model®^. 


C. Bosonic representation of charge and spin 
fermions 


Given (69), it is straightforward to rewrite (54) in terms 
of the new Bose fields 

Hb = jdx ^p, 0;)2 + 

K a 

+ E C 0 S($:/V^) 

a 

+ ^9,$:cos($:/y2) + ...|, (70) 

where — Va T 4Ka ^ — TgH 4Ka ^ ’ 


IV. LUTHER-EMERY (LE) POINT FOR SPIN 
AND CHARGE 


A particular case of the family of Hamiltonians (65) de¬ 
scribes a free theory of non-interacting gapless fermionic 
spinons and holons. This LE point for both spin and 
charge corresponds to 


= y R\{-iv'cdx -gdl-^ ... )R^ 

+Ll{iv'^dx -gdl + ... )L^ 

+Rl{—iVgdx + iCdl . )Rg 

+Ll{iv'gdx - iCdx + ■ ■ ■)Lg^ 


(71) 


Here we have included the quadratic (cubic) term in the 
holon (spinon) dispersion to emphasize the nonlinearity. 
In order to realize a Hamiltonian of this form fine-tuning 
a number of couplings is required, as can be seen by an¬ 
alyzing the stability of (71) to perturbations. 


Our spin and charge fermion fields can be bosonized 
by standard methods. Introducing chiral charge and spin 
{a = c, s) Bose fields ipa, (pa, and ignoring higher har¬ 
monics, we have 


Ra {x') 




Rql (^) 



V>*c(x) 


, (67) 


where 77 ^, ga are Klein factors fulfilling anticommuta¬ 
tion relations {ga^dp} = 2 ( 5 a ./3 = {Va.'np}, {ga,gp} = 0 . 
Bosonizing the spin and charge fermions leads to the fol¬ 
lowing expressions for the original right- and left-moving 
spinful fermions 

R^{x) oc Y[ , 

0' — C,S 

L^{x) oc n . ( 68 ) 

0' — C,S 

The new Bose fields (^* , 1 ^* are related to the usual spin 
and charge bosons (54) by a canonical transformation 


$ 


a — 



(69) 


A. Stability of the LE point 


An obvious question is to what extent the LE point is 
stable. The most important perturbations to (71) are 


Rpert 


dx 


gc.oR^RcRc^c 


+gs,i{RlLg LlRg) gg^oRlRgLgLg 
+gi{RlRc + L\L^){R\Lg -I- h.c.) 


(72) 


In addition to (72) there are other, less relevant pertur¬ 
bations. A list of the ones with scaling dimension be¬ 
low four is given in Appendix B. The ggp term in (72) 
is recognized as a mass term for spinons, and is the only 
strongly relevant perturbation. This implies that spinons 
are generically gapped, and in order to reach a LE point 
with gapless spinons fine tuning ggp = 0 is necessary. 
Assuming that this is possible, we are left with three 
perturbing operators of scaling dimension 2. While the 
5 s 0 5 c ,0 terms are scalar, the 51 term carries non-zero 

Lorentz spin. In order to assess the stability of the LE 


0 „ = v^ 0 ; . 
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point to these perturbations, we have carried out a renor¬ 
malization group analysis. In principle we need to work 
with different cut-offs for the charge and spin degrees of 
freedom. However, at one-loop logarithmic divergences 
are encountered only in the spin sector. We obtain RG 
equations of the form 


dgi 

1 

dvc 

1 

d£ 

. 9i9s,0\ 

A.'KVs 

' d£ 

2'k'^v. 

dgc.o 

1 2 

dgs.o 

:^=0 

d£ 

dl 

=- 9i, 

TTVs 

d£ 


(73) 

(74) 


where £ = ln(L/ao) and oq and L are short and long 
distance cutoffs respectively. The RG equations are easily 
integrated 


ffs.o(^) — 5s,o(^o) 
ffc.o(£) = ffc.o(£o) 


9s,0(^0) 


9i{() = gi(£o)e'‘"”=<'o) 
ffs,o(£o) 


H-io 


- 1 


(75) 


and imply the following: 


1. The spinon mass term is not produced under the 
RG flow if the bare coupling is initially set to zero. 
We have checked that this remains true at two 
loops. However, we cannot rule out that may 
be generated at higher orders and it is possible that 
setting it to zero requires fine tuning an infinite 
number of parameters in a lattice model. 

2. The coupling ^ does not flow under the RG. This 
remains true at two-loop order. Hence, to this or¬ 
der, 5s 0 needs to be fine-tuned to zero in order to 
reach the LE point. 

3. If the initial value 5s,o(£o) < 0, the coupling gi{t) 
flows to zero under the RG, while gc,o{^) flows to a 
constant value. 



Figure 2. Support in the energy/momentum plane of excita¬ 
tions with the quantum numbers of an electron or hole. For 
commensurate band fillings there is an absolute threshold that 
we take to follow the spinon/anti-spinon dispersions. Above 
the threshold the single-particle spectral function is singu¬ 
lar, and we aim to determine the threshold exponent of the 
negative-frequency part at a momentum kp — q (green circle). 


For commensurate band fillings the spectral function 
has a threshold at low energies. To be specific, we will 
consider the case < n/, in which case the threshold 
corresponds to exciting a single high-energy spinon, while 
(anti)holon excitations have vanishing energy. The corre¬ 
sponding kinematics is sketched in Fig. 2. The negative 
freqnency part of the spectral function at fixed momen¬ 
tum transfer k exhibits a threshold singularity 


A(w, k) 


0 if 0 > w > -es(fc) 

Ao|w -|- es{k)\^ if w —>■ —es(fc) . 


Here es{k) denotes the spinon dispersion. 


B. Threshold singularities in the single electron 
spectral function 


Given the low-energy Hamiltonian at the LE point 
(71), we are now in a position to derive a mobile im¬ 
parity model, valid a priori at low energies. The nsnal 
continuity arguments snggest that the restriction to low 
energies can be relaxed and the model applied to energies 
of the order of the lattice scale t. Let us focus on the mo¬ 
bile impurity model relevant for analyzing the threshold 
behaviour in the single-electron spectral function 


Aiuj, k) 


Gj-et (^, £-) 


-Im Gret(w, k), 

TT 

poo _ 

-i J dt ^ 

X (V'o|{G+i,^(£), (76) 


where |'0o) is the ground state. 


C. Threshold exponent at the LE point 

Let us focus on momentnm transfers kp — q, where we 
take 0 < q kp. Using the decomposition 

C. ~ , (78) 

we see that the relevant field theory correlator is 

pOO pOO 

H(w, kp-q)^-i dt dx 
Jo J — oo 

X (V’o|{-Rt(£, 2 ^), 7?|(0,0)}|'!/'o) 

= A<(w, kp - q) + kp - q). (79) 

Here and are respectively the positive and nega¬ 
tive frequency parts of the spectral fnnction. Using (62) 
we arrive at the following expression for the latter 

pOO pOO 

A<{uj, kp-q)^-i dt dx 

Jo J — oo 

X n -(80) 
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At the LE point we are dealing with a free fermion theory. 
Hence correlation functions of the kind required in (80) 
can be expressed as Fredholm determinants®®, but we 
do not follow this route here. Instead, we construct a 
mobile impurity model and use it to extract the threshold 
exponent. 

At the LE point spin and charge degrees are perfectly 
separated. As a consequence it is possible to construct a 
basis of energy eigenstates in the form 

W) ® \ns), (81) 

where ric^s ^-re appropriate charge and spin quantum 
numbers. The correlators required in (80) then have 
Lehmann representations of the form 

5 ]e*^-‘-*^-"|(n„|0„(0,0)|V'o)|" , a = c, s. (82) 

Ua 

The threshold singularity arises from excitations involv¬ 
ing a single high-energy spinon with momentum q plus 
low-energy excitations in the charge and spin sectors sec¬ 
tor. This means that the charge part of (80) can be cal¬ 
culated using bosonization. The bosonization identities 
(68) imply that 

(83) 

At the LE point the coupling constants Aq, Ai, A 2 in (70) 
vanish, and a simple calculation gives 

(O|(0,0)C>c(t,a;)) oc (u't - x)“5 {x^ - v'^t^)~^. (84) 

In order to work out the contribution from the spin part 
we follow Ref. 67. We decompose the spin fermions into 
low-energy and mobile impurity parts 

Rs{x) - rs[x) + , 

Ls{x) ls{x), (85) 

where xl creates a hole in the spinon band. In terms of 
momentum modes 

Rs{x) ^ I ^ e*P"R,(p), (86) 

this projection corresponds to 

rs{x) - t ^ 

J-\' 27T 

^ , A" « |g|. (87) 

Jq-A" 


Substituting the decomposition (85) into our expression 
of the Hamiltonian density (71) and dropping oscillatory 
contributions under the integral, we arrive at the follow¬ 
ing mobile impurity model 


^MiM = J dx - lidj^) 

Ot — C,S 

+ J dx xl {£s - iusdx + ...)xs ■ (88) 

Here we have introduced notations Eg = €s{kp — g) « 
<9 + + ... and ug = - ^ « < + 3Cg2. 

Next we need to work out the projection of the oper¬ 
ator Os on low-energy (rg) and impurity (xs) degrees of 
freedom. Using that rs{x), Isix) and Xs{x) are slowly 
varying fields, we can approximate the string operator as 


Osix) ^ [rg(x) [^’^)U)U(U-h.c.] 

(89) 


Here the second term is the contribution of the string 
arising from the upper boundary of integration x. By 
virtue of the presence of the strongly oscillatory factor 
^iqx expression (80) for the spectral function, the 

leading contribution to the spectral function arises from 
the part of Os{x) proportional to 

C>g(x) ... (90) 


In terms of this component we have 


A^{uj,kF — q) ^ —i / dt dx e“*((!lj(0, 0)(!lc(t, x)) 
Jo J — 00 

x(O(«)\0,0)a(La:)) . (91) 


In order to isolate the desired contribution, we expand 
the second factor in (89) 


C>(«)(x)~e *«“|rg(x) ^rl{x)xl{x) + ... -fxl(x) 1 - ^{xs(x)rg(x),rj(x)xl(x)}-f ... 


X e 


(92) 
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In order to proceed further, it is convenient to bosonize 
the low-energy degrees of freedom associated with rg, Is 
using ( 68 ): 

[bo + bida,ipl{x) -f ...] 

X (- 93 ^ 

This can be simplified further using the appropriate op¬ 
erator product expansions. In order to determine the 
threshold singularity it is sufficient to retain only the 
term with the lowest scaling dimension at low energies, 
which is 

0^^\t,x) ~ (94) 


The two-point correlator of this operator is 


( 0 (^)^( 0 , 0 ) 0 (^)(t,x)) 


2 „,, 2 , 2 x-i ,e.tSin[A"(a;-Ugt)] 


tt{x — Ugt) 


.(95) 


Sending the cutoff h" to infinity turns the last term into 
a delta function <5(a: — Uat). Substituting the resulting 
expression for (95) and the charge sector contribution 
(84) into the expression (80) for the hole spectral function 
and then carrying out the space and time integrals, we 
arrive at the following result for the threshold behaviour 


A^(bj,kF - q) o:\u! + £s\ (96) 


As expected there is a threshold singularity. The expo¬ 
nent is seen to be momentum independent. As we will 
see, this is particular to the LE point. 


V. MOBILE IMPURITY MODEL AWAY FROM 
THE LE POINT 

We now wish to generalize the above analysis to the 
Luttinger liquid phase surrounding the LE point. We 
will assume that 

1 . the spinon mass term is fine-tuned to zero, i.e. 

5 s .1 = 0 ; 

2 . the four fermion interactions in the spin and charge 
sectors are attractive, i.e. gc,o,gs,o < Oj a-nd size¬ 
able. 


= 


167r 


2K, 


+ 2Ks,{d,e: 


'H-irap — Xs (^s iUs^x) Xs i 


'biiYit xixs 


XI fo.{q)dx‘fa + fa{q)dx^* 


• (98) 


Here we have dropped all terms that do not affect the 
threshold exponent and retained the same parameteriza¬ 
tion of the impurity part of the Hamiltonian, although 
the actual values of and Ug are of course not the same 
as the LE point. The Luttinger parameter in the spin 
sector varies from Kg = 1/2 at the LE point to Kg = \ 
in the SU(2)-invariant limit. The charge Luttinger pa¬ 
rameter equals Kc = 1/2 at the LE point, and varies 
with doping and interaction strength otherwise. We note 
that close to the LE point (in the sense that go,c, 5o,s are 
small), there is an additional contribution to "Hint of the 
form xiXs cos($*/-\/2). The analysis of this case is very 
interesting (see e.g. Ref. 68 for a related problem), but 
beyond the scope of our work. The functions /a ( 5 ), /«(<?) 
as well as the parameters Va, K^, Eg, Ug depend on the 
microscopic details of the particular lattice realization of 
our field theory. We will show below how they can be 
fixed either numerically in the generic case or analyti¬ 
cally when our theory is applied to the Hubbard model. 

The mobile impurity model (98) can now be analyzed 
by standard methods®'*^’®^. The interaction between the 
impurity and the low-energy degrees of freedom can be 
removed through a unitary transformation 


If — d.^I2A'ya‘Pa{x)+'ro,Va(^)]xUx)xs{x) 


(99) 


The transformed spin impurity field equals 
dg{x) = Uxs{x)U^ 


( 100 ) 


while the chiral spin and charge Bose fields transform as 

= Uip*Jx)U'< = iflix) - 2TTXaC{x), 

‘faix) = U‘fl{x)u'< = !fl{x) + 2TrjaC{x), (101) 


where 

C{x) 
We note that 


dy sgn(a: 


y)xliy)xgiy) ■ 


( 102 ) 


Under these assumptions holons and spinons remain gap¬ 
less, and the gi term in (72) is irrelevant so that we can 
drop it at low energies. Focussing again on the single¬ 
electron spectral function, using the decomposition (85), 
and finally bosonizing the low-energy spin and charge de¬ 
grees of freedom, we arrive at a mobile impurity model 
of the form 




MIM 


= dx 


E 


"Hr, -f Himn + Td 


int 


(97) 


dxy^lix) = dx^flix) - ^T^laxl{x)Xg{x), 

dx^lix) = dx(pl,{x) + ATrjaxlix)Xgix). 


Adjusting the parameters 7 q,, 7 ^ such that 
ffA _ fu-v+ -V- \ f 7a\ 

\U) V + 


with 


2 iX„± 


2Kr, 


(103) 

(104) 


(105) 
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the impurity decouples in the new variables 




MIM 


= dx 


%!' - “ 


E njf I njf 

iLo^ ~r ^imp 
s 


,Ct — C,S 

1 


2K, 


^imp — iUgdx) dg . 


(106) 


The interaction between the mobile impurity and the 
Luttinger liquid degrees of freedom is now encoded in the 
boundary conditions of the transformed Bose fields <i)°, 
0 °, which are “twisted” by the presence of the impurity, 
see e.g. (103). The negative-frequency part of the single¬ 
electron spectral function is again given by (91), where 


Oc{x)Oi‘^\x) - .(107) 

In terms of the transformed fields this reads 
Oc{x)Oi‘>\x) - 


X e 


(108) 


The threshold behaviour of the hole spectral function can 
now be calculated in the same way as at the LE point. 


The result is 


Ag^{uj,kF - q) 


1 


|w-f eor 

where the exponent /r is given by 


(109) 


Fc,± = 


M — 1 ~ 2(r'c__|_ -I- _ -I- _|_ -I- i'g _), 

1 


2 1 


7c + 7c 


± 


2 V 2 J VW~c^ 


Ic-lc- 


1 1 

71 


1 


± 


: [is - Is] ■ ( 110 ) 


■'.,± = VTr'+^' + 2V2J^vw; 

As 7 a, 7a are functions of q, the threshold exponent 
is now generally momentum dependent. However, it is 
shown in Appendix C that spin rotational SU(2) symme¬ 
try in the limit Kg —>■ 1 enforces the particular values 


1 


Is = Is = 


4^2’ 


( 111 ) 


for any value of q. 


A. Relation of 7^, 7^ to finite-size energy spectra 

An obvious question is whether there is a way of di¬ 
rectly determining the parameters 7 a, 7a for a given mi¬ 
croscopic lattice model. To that end, let us consider the 
spectrum of our mobile impurity model on a large, finite 
ring of circumference L. The mode expansions of the 
Bose fields (^* , are 


n—1 






ct,R,n 




( 112 ) 


Here Q*, Q* , (/?* q and (/?* q are zero mode operators with 
commutation relations 


Ko,Q:] = -47r* = -Ko,o:]- 


(113) 


The eigenvalues qa, qa of the zero mode operators Q*, 
Q* depend on the boundary conditions on the fields (/?*, 
(^a, which on general grounds will depend on whether 
or not a mobile impurity is present. In presence of the 
impurity the finite-size spectrum of (98) has the following 
structure 

(114) 


Here Ei^p is the contribution of the impurity to the en¬ 
ergy. On general grounds it will have the following ex¬ 
pansion in terms of system size 


= eZ + iE\^^ + o(L-i). 


p(i) 


(115) 


The other contribution to (114) arises from the Luttinger- 
liquid part of the theory. Applying the mode expansions 
to the transformed Hamiltonian (106) we obtain 


E — Eqs + Aimp -|- AAll + o{L ^). 
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ASll = E 


2TrVn 




Qa “t" Qa 
in 


- 7a + 7c 




Qa Qa 
in 


la - Ic 




M- 


n>0 


• (116) 


Here Qa, Qa ^-i^d are “quantum numbers” charac¬ 

terizing a particular low-energy excitation. Their quan¬ 
tization conditions depend on the boundary conditions 
for the spin and charge Bose fields (112) in presence of a 
high-energy mobile impurity. These can be worked out 
by considering the “minimal” excitation that can be made 
in the sector where the mobile impurity is present. This 
sector is reached by acting with the operator 

0,{x)0\^\x) = (117) 


on the ground state |^/’o)- The latter is characterized by 
being annihilated by aa,R,n, <ia,L,n, Qa Qa- The 
quantum numbers of the “minimal” excitation are found 
by noting that 

Q*aO,{x)Oi‘^H^)\^o) = [Q*a,0,{x)Oi’^Hx)]\^Po) 

= Q^°'>0,{x)Oi<^H^Mo), (118) 


where 


o(o) = ^ do) ^ JL „(o) = _do) ^ _!L 

Qc ^.Qc ^,Qs Qs ^2-^ ^ 


Let us denote the lowest energy state with these quantum 
numbers by 






( 120 ) 


Higher excited states in the “impurity sector” can be ob¬ 
tained for example by making particle-hole excitations 
on top of the state (120). The energies of such states are 
given by (116) by choosing (119) and in addition taking 
some of the M)^a different from zero. Crucially, the val¬ 
ues of 7 q,, 7 o! are the same as for the state (120). For prac¬ 
tical purposes, states with the same momentum as ( 120 ) 
but different spin and charge quantum numbers might 
be of particular interest as they are lowest energy states 
in certain sectors of quantum numbers and can therefore 
be more easily targeted in DMRG computations. The 
quantum numbers for such states are 



where TOq, rha are integers. This follows from the mode 
expansion and the requirement, that the bosonized ex¬ 
pressions for R^{x), Ri{x), Li^{x), Li{x) must be single 
valued. Using that at low energies the charge and spin 
densities and currents are given by 

Pa{x) = ) 

Try 8 

ja{x) = -(x) , (122) 

Try 8 


• (ttt-c -I- mc)l2 is the difference in charge (parti¬ 
cle number) between the excitation and the state 
( 120 ); 

• {rjic—fric )/2 is the number of charge fermions trans¬ 
ferred from the right-moving to the left-moving 
branch; 

• (ms+fhs)/2 is the difference in the number of down 
spins between the excitation and the state ( 120 ); 

• {rris — rha) 12 is the number of spin fermions trans¬ 
ferred from the right-moving to the left-moving 
branch. 

The values of 7 ^, 7 q, can now in principle be extracted 
by numerically computing finite-size energy levels by a 
method such as momentum-space The pro¬ 

cedure is outlined in the following. 

(1) The velocities Va and Luttinger parameters Ka are 
bulk properties and can be determined by standard meth¬ 
ods. We will assume them to be known quantities in the 
following. We further denote the number of particles and 
down spins in the ground state by IVqs and Mgs respec¬ 
tively. 

(2) We then numerically compute the lowest excitation 
with momentum q and quantum numbers N = Nqs ~ Ij 
M = Mgs ~ 1- Its energy is 

+ AF;LL(gf + o(L-i).(123) 

(3) Next we compute the lowest excitations with momen¬ 
tum q but different values of N and M. For example, 
choosing 


TV = TVgs - 1 , M = Mgs - 2, 

gives the excited state characterized by 

rric = —fhc = 1 , rris = —2 , fhs = 0 . 

The zero-mode eigenvalues follow from (121) 


Stt 

In 


72 ’ 72 


71) = -— 71)=- — 

72 ^ ^ 72' 

The corresponding energy is 

E-1-2= EqS + -Uimp + 


(124) 


(125) 


(126) 


+AEi,l^{qi^\qP,qi^\qW) + 0(^-7 (127) 


we can identify 
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Here we have asserted that the change in the impurity 
contribution to the energy is of higher order in L~^. This 
has been shown for the case of the Hubbard model using 
methods of integrability in Ref. 53. We believe that 
this continues to hold true in general, because is 
sensitive only to the values of the parameters 7 a, 7 a, 
which are the same for all excitations we are considering. 

The point is that by considering energy differences like 

E_i_2 - 

(128) 

we obtain a set of equations in which the only unknown 
parameters are the 7a, 7a- This provides a numerical 


method for determining them in a general lattice model. 
For integrable theories like the Hubbard model, analyti¬ 
cal techniques are available and we discuss this case next. 


B. Hubbard model 

The finite-size spectrum in presence of a high-energy 
spinon excitation was calculated for the case of the Hub¬ 
bard model using the Bethe Ansatz solution in Ref. 53. 
The result for the lowest excited state above the spinon 
threshold is 


if = ifGS - 


-k 


27ru„ 


(AW - iV™P)2 
8 W 


2W D, - + 


D, 


2ttVs 

^17 


(AW -jAW -i) 

2 


o(L-i),(129) 


where 


w = w = 0 , AW = -1 , AW = 0. (130) 


The contribution —es{A^) — ^e(,(A^)i5A^ is the finite- 
size energy of the impurity. The velocities Va, Ka as 
well as the quantities and ii™P are expressed in 

terms of solutions to coupled linear integral equations, 
and in practice are easily calculated numerically with 
very high precision. By construction the quantum num¬ 
bers (130) correspond to our minimal excited state (120). 
By matching the Luttinger liquid part of the energy to 
(114) we then obtain the following results for the param¬ 
eters 7 a, W 


7c + 7c 


1 

2^2 


Is =78 = 


- v^W“P, 
1 

4v^' 


lc-lc = - 


N'nnp 


(131) 


One subtlety to keep in mind when making contact be¬ 
tween the Bethe Ansatz calculation and (114) is that the 
former refers only to highest weight states of the SU(2)(g) 
SU(2) symmetry algebra of the Hubbard modeR^A2 p)g_ 
scendant states need to be taken into account separately. 
Substituting (131) in the expressions for the quantities 
r'a (110), we find 


1 -I- 

= 0, T . (132) 

Finally, the threshold exponent is obtained from (110) 

M = 1 - (133) 


This agrees with what was found in Ref. 53 using the ap¬ 
proach of Schmidt, Imambekov and Glazman®^’^^, as well 
as with the exponents reported previously in Refs. 31-33. 


1. Other excited states 

The Bethe Ansatz result (129) can be applied to other 
excited states as well. In particular, the holon plus three 
spinon excitation considered in Ref. 53 gives rise to an 
excitation with the same momentum and quantum num¬ 
bers 


AN, = -1 = AN, , W = -1 , W = ^- (134) 

The corresponding state in our mobile impurity model 
has quantum numbers (125), and the energies calculated 
from (114) and (129) agree as they must. 


VI. RELATION TO THE APPROACH OF 
SCHMIDT, IMAMBEKOV AND GLAZMAN 

The method of Refs. 51 and 52 is based on a dif¬ 
ferent prescription for defining fermionic quasiparticles 
in the charge and spin sectors. We now summarize the 
main steps of this approach and compare them to our 
framework. The starting point is the standard bosonized 
description (53) of the spinful fermion model under con¬ 
sideration. One then introduces the chiral components 
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4^a t 0a by 

$„(a;) = + <^Q,), (135) 

Qa{x) = ^={(j)a- ^a)- (136) 

V -^a 

These chiral bosons diagonalize the spinful Luttinger 
model (53) in the form 

y-Lhix) = X! ^ [{dx4>af + [dx^af] ■ (137) 

OL — C,S 

Next, in analogy with the spinless case in Eqs. (41) and 
(42), one defines the quasiparticle operators 

Rc{x) ^ La(a;) ~ (138) 

As in the spinless case, this prescription removes the 
marginal interactions between quasiparticles, rendering 
them asymptotically free in the low-energy limit. How¬ 
ever, these quasiparticles cannot be identified with the 
finite energy elementary excitations in integrable mod¬ 
els, because they carry fractional quantum numbers: 

[qa,,RUx)] = V^Rl [qo.,LUx)] = V^Ll 

(139) 

The only exception is the LE point = 1/2, where 
the operators Ra , La in fact carry the same quantum 
numbers as our spin and charge fermions. 

Despite the lack of correspondence with long-lived ex¬ 
citations in (nearly) integrable models, the approach of 
Refs. 51 and 52 allows one to compute threshold expo¬ 
nents, as long as the parameters in the effective impurity 
model are adjusted appropriately. The situation is analo¬ 
gous to the two ways of determining threshold exponents 
discussed above in the spinless fermion case. 

In order to facilitate a direct comparison with our ap¬ 
proach, we briefly review how to express the electron op¬ 
erator in order to calculate the threshold exponent in the 
single-particle spectral function within the approach of 
Refs. 51 and 52. We assume again that the lower thresh¬ 
old of the support corresponds to an excitation with a 
single high-energy spinon. First, the right-moving spin 
quasiparticle in Eq. (138) is projected into low-energy 
and impurity subbands: 

Rs{x)r^fs{x) + e~''’^'^x\{x). (140) 

One then rewrites the electron operator in terms of the 
quasiparticles defined in (138). As this differs from ours, 
cf Eq. (67), the string-operator part in the expression for 
the electron operator is also different 

Oi—C^S 

X (141) 

The next step is to write down an effective impurity 
model, analogous to Eq. (98), but using as the im¬ 
purity. After performing a unitary transformation that 


removes the coupling between Xa and the low-energy 
modes, the electron operator becomes 

R^{x) - 

X g-1 [ (7^ -47I) 7^ -47') 0c] 

xe-i[(v^-7fc-470^^], (142) 

where dg = UxsU^ is the free impurity field for the quasi¬ 
particle with fractional charge, and 7 /,, 7 ), are the param¬ 
eters of the unitary transformation, which are not the 
same as 7 a, 7 a discussed in Section V. 

At the LE point, we set Ka = 1/2 and 7 /^ = 7 ),, = 0 
and Eq. (142) reduces to 

R^{x) ~ (143) 

This result agrees with the refermionization in Section 
IVC, since at the LE point da = da = Xa^ 4>a = ‘Pa, and 
4'a = Pa', thus, the two approaches coincide. 

Moving away from the LE point, the expressions for 
physical operators in terms of impurity and low-energy 
fields will in general be different in the two approaches. 
However, the results for the edge exponents are still con¬ 
sistent because the difference in string operators can be 
accommodated by the parameters of the unitary trans¬ 
formation and by imposing proper boundary conditions 
on the bosonic fields. For instance, imposing SU(2) sym¬ 
metry in the approach of Refs. 51 and 52 leads to the 
requirements 

7 ( = V2-1, 7 / = 0 . (144) 

This should be contrasted with Eq. (111). 

VII. REALIZING THE LE POINT 

We have seen that a good starting point for under¬ 
standing threshold singularities in dynamical response 
functions is the LE point for both charge and spin. In 
section IV we considered properties of the LE point in the 
field theory limit. An obvious question raised by these 
considerations is whether it is possible to realize the LE 
point in practice in a lattice model of interacting spinful 
fermions. We now investigate this issue in some detail 
and present a number of preliminary results. 


A. Lattice model 

As discussed in section IV A, realizing the LE point at 
sufficiently low energies in an extended Hubbard model 
of the kind (47) requires the fine-tuning of (at least) four 
parameters: 

K, = Ka = ^, ga.i = 51 = 0. (145) 
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For the Hubbard model in zero magnetic field, spin ro¬ 
tational symmetry fixes = 1, while Kc varies with 
both band filling and interaction strength U. In partic¬ 
ular, it is well known^^’^^ that Kc = | is obtained in 
the U ^ oo limit of the Hubbard modeb^’^^. Values 
Ks < 1 can be realized by adding spin-dependent inter¬ 
actions that break the spin SU(2) symmetry, but retain 
spin inversion symmetry. The latter is crucial for avoid¬ 
ing marginal interactions between spin and charge sectors 
that lead to a more complicated conformal spectrum in¬ 
volving a dressed charge matrix^’^^’^^. 

A minimal lattice model that may allow us to fulfil the 
conditions in (145) is 

L-l L 

j=^l (7 j^l 

2 L-r L-l 

+E E E ■ ( 146 ) 

r=l j = l j=l 

In anticipation of the DMRG computations of energy lev¬ 
els reported below we have imposed open boundary con¬ 
ditions. In the following we set t = 1, i.e. measure all 
energies in units of the hopping parameter. The idea 
is then to try to adjust the four interaction strengths 
U, Vi ,2 and Jf in such a way that (145) are achieved. 
In order to ascertain the low-energy properties of (146), 
we compute the energies of the ground state and sev¬ 
eral low-lying excited states for a quarter-filled band, and 
compare the results to expectations based on Luttinger 
liquid theory. We choose to work at quarter filling in or¬ 
der to simplify finite-size scaling analyses. As alluded to 
in Appendix A, working at commensurate fillings induces 
additional Umklapp interactions. In the case at hand this 
corresponds to the presence of an additional perturbation 
f dx cos(24>c). However this term has scaling dimension 
SATc and is therefore highly irrelevant for Kc >1/2. We 
therefore discard it in the following analysis. 

Some insight into how the parameters Kc,Ks and i 
depend on U, Vi ^2 and Jf can be gained by bosonizing 
the interactions at weak coupling. Using (48) we obtain 
in leading order 


rijUj+i 

njnj +2 


QZ QZ 

W W-si 


A 

4 1 


of (x) + 20f (a;)-Of (x) 

-20f (x)(x) , (147) 
of (^) + 20f (x)l , (148) 


A 

A 


2 L 


30f ^ (x) -k 60f (x) + 0 ^ 2 '’ (x) 
-k20f (x) - 87/tr?;ry;%of (x) , (149) 


of (x) + 20f (x) 


(150) 


Here x = jop, A is a dimensionful amplitude, and 

of = cos , 

O^^^ = {dx(psf + {dx'^sf , 

O^^^ = dx<fsdx<fs , 

= {dx'^cf + {dx^cf , 

of = dx^Psdx^s- (151) 

At weak coupling the spin-charge separated Luttinger liq¬ 
uid at low energies is therefore perturbed by 



5 

dx of (x) = 

3 = ^^ 


(152) 


where the Vj are proportional to linear combinations of 
U, Vi, V 2 and Jf. The effects of dH 2 and 6 H 4 are to 
renormalize the spin and charge velocities respectively, 
while SH^ and 6 H^ change the values of the Luttinger 
parameters K^ and Kc- We must pay special attention 
to the perturbation SHi'. as discussed in Appendix B, 
this operator is marginal at weak coupling, but gives rise 
to the relevant spinon mass term as we approach the LE 
point. Our objective is to adjust U, Vi, V 2 and Jf in 
such a way that iFc,s are reduced towards 1/2, while the 
coupling vi of the spinon mass term remains very small. 
The bosonization results (150) suggest the following pre¬ 
scription for achieving this at weak coupling: 

• Increase V 2 in order to make \i'i\ very small; 

• Then adjust Vi and Jf in order to drive ATc.s to¬ 
wards 1/2. Importantly, the bosonization results 
(150) indicate that at least at weak coupling this 
does not produce sizeable contributions to 5Hi. 

Our analysis below is guided by these considerations, 
even though we are not operating in the weak coupling 
regime, in which (150) are applicable. 


B. Finite-size spectrum of unperturbed Luttinger 
liquid with open boundary conditions 


A standard procedure for determining the Luttinger 
parameters Kc and Kg is to compare the finite-size spec¬ 
trum predicted by Luttinger liquid theory with the low- 
energy spectrum calculated numerically for a given lattice 
model. The finite-size spectrum relative to the ground 
state of a Luttinger liquid with open boundaries is given 

by 


A£;(5",AW,{m/,m|}) 


7rUs(5^)^ TTVciANc)'^ 
KgL ^ 4:KcL 


TTVcdANc 

2 KcL 


00 

+ E E (^^^) 

Oc — C,S £—1 


Here and AiV^ are quantum numbers in the spin and 
charge sectors, respectively, associated with the global 













17 


U(1)0U(1) symmetry. The value of measures the 
change in the total magnetization and ANc = N — Nq 
measures the change in the total nnmber of electrons with 
respect to a singlet ground state with Nq electrons. As 
nsnal, the values of and ANc are constrained by the 
selection rnle that 2S^ must be even (odd) if ANc is 
even (odd). The parameter d in Eq. (153) is a dimen¬ 
sionless constant that depends on the definition of the 
chemical potential to order 1/L. The parameters and 
ml are non-negative integers that count the nnmber of 
low-energy particle-hole pairs in each sector. We note 
that for open bonndary conditions the total momentnm 
is not conserved and there are no qnantnm nnmbers as¬ 
sociated with current excitations (i.e. with transferring 
particles between Fermi points). Thns (153) should not 
be confused with the spectrum in (116), which is valid 
for periodic bonndary conditions. 

The spin and charge velocities Vs and Vc can be ex¬ 
tracted by matching the energies of the lowest excitations 
with = ANc = 0. Assuming Vs < Vc (which can be 
verified by analyzing excitations with different quantnm 
numbers), the first excited state corresponds to mf = 1 
and has energy 

-TTII 

AE{0,0,{m^f> = = 0}) =—j^. (154) 

Lj 

Likewise, if Vg < Vc < 2vs, the second excited state cor¬ 
responds to ml = 1 and has energy 

TTIJ 

AE{0,0,{mt = 0,m1 = 5i^i}) = (155) 

Jj 

Having determined the velocities, one can obtain the 
Luttinger parameters by analyzing low-lying excitations 
that change the quantum numbers and ANc- We 
adopt the short-hand notations 

AEoiS^ANc) = AE{S^,ANc,{m'i = = 0 Vi}) 

(156) 

for the lowest energies in each sector of fixed and 
ANc- To isolate the dependence on Kc, we consider ex¬ 
citations with ANc = i2 and = 0. The dependence 
on the unknown constant d can be eliminated by taking 
the combination 

^[AAo(0,2) + AAo(0, -2)] (157) 

where we recognize k as the compressibility of the Lut¬ 
tinger liqnid^^. Analogonsly, Kg can be determined nsing 
the relation for the finite-size spin gap 

§AAo(1,0) = ^. (158) 

The right-hand side of Eq. (158) is equal to the inverse 
spin snsceptibility. 

The procedure described above is standard for Lut¬ 
tinger liquids that are only perturbed by (strongly) ir¬ 
relevant operators. In our case, however, we must con¬ 
sider the effects of the relevant operator cos<i>s, which 


generates the spinon mass, as well as the interaction 
9a;‘he cos 4)s, which is marginal at the LE point and only 
weakly irrelevant in its vicinity. These perturbations 
introduce corrections to (154), (155), (157) and (158), 
which affect the finite-size scaling analysis. Since cos4>s 
is the leading pertnrbation, we first foens onr efforts on 
fine tuning it to zero, as we disenss in the next subsection. 


C. Fine tuning the spinon mass 

As anticipated in Sec. VII A, onr strategy for fine tun¬ 
ing the spinon mass to zero starts by snppressing the 
conpling constant of the pertnrbation SHi. According 
to (149), for fixed U > 0 this can be done efficiently by 
increasing the next-nearest-neighbor interaction V 2 > 0 . 
In this process, we keep Vi = Jf = 0, so the SU(2) 
symmetry is preserved. As we keep increasing V 2 , the 
marginal coupling constant will change sign at some 
critical value Vf. Beyond this point the perturbation 
becomes marginally relevant and the system nndergoes 
a Berezinskii-Kosterlitz-Thonless (BKT) transition to a 
spin-gapped phase, analogons to the dimerization transi¬ 
tion in the J 1 -J 2 spin chain^^. 

Becanse of the exponentially small value of the gap in 
the vicinity of the critical point it is difhcnlt to pinpoint 
V 2 by means of a gap scaling analysis. A better approach 
is to determine the critical point by searching for a level 
crossing in the spin excitation spectrum in the nentral 
sector = ANc = 0^®. The idea is that in a Lnttinger 
liquid the state with m| = 1 should be degenerate with 
the state mf = 2 ; however, in practice the degeneracy is 
lifted at order (LlnL)“^ dne to the marginal perturba¬ 
tion. Exactly at the critical point, the conpling constant 
for cos(4)s) vanishes, and V 2 can be identified from the 
level crossing between m^ = 1 and mf = 2 states. 

D. Analysis of finite-size excitation energies 

After tuning the conpling constant of SHi to zero, we 
proceed to adjnsting Vi, and Jf in order to achieve valnes 
Kc,g « 5 . This has to be done while making snre that 
the coupling constant in (65), or eqnivalently Ai in 
(54), is held at zero as we approach the LE point. To 
determine the values of Kc,g, and also to ascertain that 
we remain in the Luttinger liquid phase as we vary the 
parameters of the lattice model, we resort to an analysis 
of the finite size energy levels of the ground state and 
low-lying excitations. For this purpose, we now refine 
the expressions presented in Sec. VIIB by inclnding the 
effects of the leading perturbations in the vicinity of the 
LE point. 

An inherent difhcnlty enconntered when approaching 
the LE point is that corrections to the Luttinger liqnid 
form (153) of the finite-size energy spectrum become in¬ 
creasingly complicated. To see this let us consider the 
bosonized form (54) of our Hamiltonian. In the regime 
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i < Ks^c < 1 we expect the structure of the energy dif- similarly. For example, for Vs < Vc, the size dependence 
ference AEq{S^, 0 ) to be of the form of the first excited state in the neutral sector is given by 


Aifo(5^0) 


ai TTVsjS^)'^ bi 

KsL 

-I-1-1-r- 

f^AKs-l j^6K^-2 ■ ■ ■ ’ 


(159) 


where we have defined L = L + 1. The origin of the 
various contributions is as follows: 


• The oi term arises from first order perturbation 
theory in the spinon mass term Ai f dx cos<i>s. As 
a consequence of our fine-tuning we expect Oi to be 
quite small, so that higher orders of perturbation 
theory can be neglected. 


Aif(0,0, {ml = = 0}) 


a? TTVs M 

— - “h —- 

L 


L^K.-l 


+ ... (162) 


By computing the energy difference (162) numerically for 
a range of system sizes, using of and b\ 2 as fit param¬ 
eters, and fixing Kg to be the value obtained from the 
analysis of 5{L) in Eq. (161), we may extract the spin 
velocity Vg- 


3. Determining the charge velocity Vc 


• The bi term arises from first-order perturbation 
theory in A 2 / dx cos <l>s (which generally is 
non-zero for open boundary conditions). 

• The 62 and 63 terms arise from second- and third- 
order perturbation theory in A 2 / dx cos <i>s re¬ 
spectively. Here we need to consider higher orders 
in perturbation theory because the bare A 2 has not 
been fine tuned and is not guaranteed to be small. 

We see that the size dependence becomes increasingly 
complex as Kg tends towards 1 / 2 , which complicates the 
analysis of energy levels. A good way of achieving an ac¬ 
curate description of finite-size energies would be through 
(two-loop) renormalization-group-improved perturbation 
theory around both the LE point and the weak-coupling 
limits. However, as this is quite involved, we content 
ourselves with the simpler form (159) (and stress again 
that our numerical analysis is to be considered as prelim¬ 
inary) . 


1. Determining the Luttinger parameter Kg 


Equation (159) provides us with an alternative way 
of determining Kg by considering the system-size depen¬ 
dence of the quantity 

^(L) = iAEo(2,0)-AEo(l,0). (160) 

Using the expression (159) we obtain 


5[L) 


ai 


+ 


jr2if,-l j^AKg-1 


+ .... (161) 


By fitting 5{L) computed numerically for a range of sys¬ 
tem sizes to (161), we obtain a value for Kg. 


2. Determining the spin velocity Vg 

Given Kg, we may determine Vg from the size depen¬ 
dence in (159). Alternatively we can consider other ex¬ 
cited states, whose finite-size energies can be analyzed 


The structure of finite-size corrections to energy levels 
of charge excitations is somewhat simpler because some 
of the terms in first-order perturbation theory vanish. 
For Vg < Vc < 2vg, the energy of the second excited state 
in the neutral sector scales as 


AE{0, 0, {ml = 0, m/ = J^p}) 

TTVr bn 

= -A + 2-^ ^ 

L L^Kg-i 


(163) 


Using the results for Kg from the analysis described in 
VHD 1, Eq. (163) provides us with a means of determin¬ 
ing Vc through a finite-size scaling analysis. 


Determining Kg 

Finally, the energy difference related to the compress¬ 
ibility (157) is found to have the following size depen¬ 
dence 

AEo(0,2) + AEo(0,-2) ttvc 6 ^ 

16 8 KcL 

= fiL). (164) 

We can use (164) together with Kg and Vc determined in 
sections VH D 1 and VH D 3 respectively, to fix the Lut¬ 
tinger parameter Kc through a finite-size scaling analysis 
of f{L), taking b^ to be a fit parameter. 


E. Numerical results 

We now turn to the numerical implementation of the 
method set out in the previous subsection. We performed 
DMRG^®’™ computations on lattices with up to L = 216 
sites, keeping up to 3000 states and running up to 36 
finite-size sweeps. In case of the Hubbard model we have 
checked the energies obtained in this way against exact 
Bethe Ansatz results and found the relative errors to be 
of order 10 “®. 

We work with model (146) at fixed U = 3 and search 
for the LE point by varying Vi, V 2 and Jf. The first step 
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is to increase V 2 keeping Vi = Jf = 0. In Fig. 3 we show 
the first four excitations in the neutral sector for L = 64 
as a function of V 2 . By tracking the evolution of the 
energies starting from V 2 = 0, i.e. the Hubbard model, 
we are able to identify the quantum numbers 
for each of these states. The two nearly degenerate spin 
descendant states of interest correspond to the third and 
fourth excited states. From the crossing of these two 
energy levels we estimate V 2 ~ 1 . 1 . 
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Figure 3. Energies of the four lowest excited states in the 
neutral sector = ANc — 0 for U — 3, Ji = Vi — 0 and 
L = 64 as a function of F 2 . The critical value Vfy is identified 
through the crossing of the two nearly degenerate states in 
the (2s, 0) multiplet. Lines are guide to the eye. 

Next, we vary Vi and Jf to bring Kc and Kg close to 
1/2. After searching in parameter space, we settle for the 
particular parameter set 

[/ = 3 , Vi= 0.85 , 14 = 1.1 , Ji = 5.2. (165) 

The analysis presented below suggests that this corre¬ 
sponds to Luttinger parameter values of Kg = 0.655 and 
Kc = 0.500. This is probably as close to the LE point as 
one can get without a more precise theoretical description 
of finite-size energies based on a renormalization-group- 
improved perturbation theory analysis. 

Figure 4 shows numerical results for the quantity 6 {L) 
defined in Eq. (160). The data have been fitted using 
(161). The resulting best-fit estimates are 

Ks = 0.655, (166) 

di = 0.009, 61 = 0.493 and 62 = —5.131. The quality of 
the fit is visibly excellent (the residuals are ^ 10“®). The 
numerical value for di is very small, confirming that we 
have almost succeeded with fine-tuning the spinon mass 
term to zero (on the scale set by the system sizes we 
consider). 

In the next step we determine the spin velocity using 
(159) and retaining the bi and 62 terms. Setting Kg = 



Figure 4. Crosses: numerical results for S{L) as defined in 
Eq. (160). The solid line is the best fit to Eq. (161) and is 
obtained for Ks = 0.655, ai = 0.009, 61 = 0.493 and &2 = 
-5.131. 


0.655 we obtain 


Vg « 1.707. (167) 

This value is consistent with the result obtained by con¬ 
sidering AE{0,0,{m^ = = 0}) in Eq. (162). 

We have verified that the energy level corresponds to 
the first spin descendant state by tracking the tower of 
lowest-lying energies in the neutral sector along a path 
in parameter space connecting the Hubbard model to the 
point (165). 

Having determined Vg and Kg, we now turn to the 
charge sector. In Fig. 5 we present numerical results for 
the quantity f{L) defined in Eq. (164), together with 
a fit to the functional form posited in (164). The spin 
Luttinger parameter is fixed as Kg = 0.655. The best-fit 
estimates are 

171 ) 

and 63 = 0.271. The residuals are of order 0(10“®). 

Finally, in Fig. 6 we present numerical results for 
AE(0,0, {mj = 0, to/ = < 5 ^, 1 }). Fitting the data to (163) 
with Kg = 0.655 results in estimates 

Vc « 2.543, (169) 

and 6 / = —1.191. Combining (168) and (169) we con¬ 
clude that 

Kc « 0.500. (170) 


F. Priedel oscillations 

The analysis of finite-size energy levels presented above 
is clearly rather involved, and an independent check on 
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Figure 5. Crosses: numerical results for f{L) as defined in Figure 7. Crosses: local density for a quarter filled band and 

(164). The solid line is a best fit to (164) with Ks = 0.655. U = 3, Vi = 0.85, V 2 = 1.1, Jf = 5.2. The solid line is a fit 

to eqn (171) with Ks = 0.65 and Kc = 0.5. 



In Fig. 7 we compare the prediction (171) to DMRG 
results for a quarter filled band and system size L = 
208. We fix the values of the Luttinger parameters to 
Ks = 0.65 and Kc = 0.5 and use the amplitudes dj and 
phase shifts pj as fit parameters. The agreement is very 
good except near the boundaries. The best fit is ob¬ 
tained for di = —0.05, ^2 = 0.83, d^ = 0.01, pi = —7r/4 
and p 2 = —7r/2. The good agreement between the ex¬ 
pected behaviour (171) and the DMRG results provides 
a consistency check on the values of iCc,s extracted from 
the analysis of finite-size energy levels. 


VIII. SUMMARY AND CONCLUSIONS 


Figure 6. Crosses: numerical results for AE{Q,Q,{ml = 
0,= 5 ^, 1 }). The solid line is a best fit to (163) for 
Ks = 0.655.’ 


the results for Kc^s would clearly be very useful. Such a 
check is provided by analyzing Friedel oscillations of the 
charge density on a chain with open boundary conditions, 
cf Ref. 77. We summarize the main steps of how to calcu¬ 
late the charge density for open boundary conditions in 
(perturbed) Luttinger liquid theory in Appendix D. For 
a quarter filled band we obtain 


1 


di 




2 " {Dj)(.k.+Ks)/2 

cos([^+|fj]j+2p2) , 

+d^- - — -■ + (171) 




where Dj denotes 


2{L -b 1) . , TTJ N 

Dj = - - sin (--- 

^ TT ^L+ 1 ' 


(172) 


In this work we have developed a new approach to 
deriving mobile impurity models for studying dynam¬ 
ical correlations in gapless models of spinful fermions. 
Our construction is based on the principle that our mo¬ 
bile impurity must carry the quantum numbers of a 
holon/antiholon (charge ±e, spin 0) or a spinon (charge 0 
and spin ±1/2). For the case of integrable models like the 
Hubbard chain, it is known from the exact solution that 
elementary excitations at finite energies are (anti)holons 
and spinons with precisely these quantum numbers. In 
integrable models these excitations are stable (i.e. do 
not decay). Breaking integrability is expected to render 
their lifetimes finite, but leaving their spin and charge 
quantum numbers intact. 

To facilitate our construction, we first derived a repre¬ 
sentation of spinful nonlinear Luttinger liquids in terms 
of strongly interacting fermionic holons and spinons. At 
a particular Luther-Emery point for spin and charge, 
holons and spinons become noninteracting. Using this 
as our point of reference, we derived a mobile impurity 
model appropriate for the description of threshold singu- 
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larities in the single-particle spectral function for a gen¬ 
eral class of extended Hubbard models in their Luttinger 
liquid phase. 

Our construction differs in important aspects from pre¬ 
vious work by Schmidt, Imambekov and Glazman^^’®^. 
However, we demonstrated explicitly how and why re¬ 
sults for threshold exponents obtained in the two ap¬ 
proaches coincide. 

Finally, we presented a preliminary analysis of the 
question of how to realize, in the low-energy regime, the 
Luther-Emery point for spin and charge in a lattice model 
of interacting spinful fermions. We showed that the struc¬ 
ture of allowed perturbations to the Luther-Emery point 
is such that fine tunings of various interactions is re¬ 
quired. Achieving these fine tunings is very delicate, and 
we discussed in some detail what problems one encoun¬ 
ters. 

Our work raises a number of interesting questions that 
deserve further attention. First and foremost, further nu¬ 
merical studies are required in order to identify a param¬ 
eter regime in an appropriate extended Hubbard model 
that, at least approximately, realizes the Luther-Emery 
point. Second, it would be very interesting to imple¬ 
ment the numerical procedure we proposed for determin¬ 
ing the parameters 70 ,, 7 ^ characterizing the threshold 
exponents. One first might want to reproduce the known 
exact results for the Hubbard model, before moving on 
to non-integrable cases. Third, in close proximity to 
the Luther-Emery point the mobile impurity model in¬ 
volves an additional marginal interaction that cannot be 
removed by a unitary transformation. It would be inter¬ 
esting to analyze its effects on the threshold exponents. 
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1 . Symmetries 

Our lattice models of interest are invariant under var¬ 
ious symmetry operations. These symmetries are inher¬ 
ited by the bosonic low-energy description (54) and we 
now discuss their realizations. 


a. Spin flip symmetry 


The lattice models of interest are invariant under ex¬ 
change of up and down spins 


0,T ^ Od- 

is realized at the level of the bosonic fields as 




(Al) 


(A2) 


b. Translational Invariance: 


The translation operator acts as 

T^Cj^aT = Cj + i^a- (A3) 


We can implement this on the level of bosonic fields by 
imposing the transformation properties 


/ <(Sc(a;)\ _ ^ f Tc{x) - 2kFao\ 

\^<^c(a;) j \Tc{x) - 2kFao) 

ff 


(A4) 


We note that this transformation works for all higher har¬ 
monics in (48). Translational invariance of the Hamilto¬ 
nian then implies that it generically may not contain any 
vertex operators of the charge boson. Exceptions to this 
rule occur at commensurate fillings 


p 

kpOQ = TT- , p, g S N. 

(A5) 

q 


Here operators of the form 


cos(|$c), sin(|$c) 

(A 6 ) 


are allowed to occur. As long as q > 4 such operators are 
highly irrelevant and do not play a role in the following 
discussion. 


Appendix A: Irrelevant perturbations to the 
Luttinger liquid Hamiltonian 

One way of working out the allowed irrelevant per¬ 
turbations to the Luttinger liquid Hamiltonian is by us¬ 
ing symmetry considerations. Our starting point are ex¬ 
tended Hubbard models of the kind (47). For the pur¬ 
poses of this appendix, we will assume all interactions to 
be small. 


e. t/(l) (8) [/(I) Invariance 

By this mean we that the Hamiltonian commutes with 
particle number and the z-component of total spin 

[H, S^] = \H, iV] = 0. (A7) 

This symmetry implies that no vertex operators involv¬ 
ing the dual fields 0 ^, Og are allowed to occur in the 
expression for TL. 


22 


d. Site Parity 

The reflection symmetry acts on the lattice fermion 
operators like 

Pcj^aP = (A8) 

We see that P can be realized in the fleld theory as 
Pipa{x)P = -(pa{-x) , PLpa{x)P = -ipa{-x) , 

PriaP = fja- (A9) 

Crucially, parity acts on the unphysical Klein degrees 
of freedom as well. Moreover, we obtain the following 
constraints on the amplitudes An^m in (48) 

An^m — Ai — n, — l — m- (-^40) 

2. Dimension two operators allowed by symmetry 

We now list all symmetry-allowed perturbations to the 
Luttinger liquid Hamiltonian for noninteracting spinful 
fermions with scaling dimensions 2, 3 and 4. Such per¬ 
turbations will be of the form 

J dx 0{x), (All) 

3 

where Xj are coupling constants and Tj are products of 
Klein factors. The only non-trivial combination of Klein 
factors allowed to appear is in fact 

(ai 2 ) 

This is because the all terms in our Hamiltonian (47) are 
of the form 

• (-^43) 

Expressing the lattice fermion operators in terms of Bose 
fields by (48) , and then imposing that for a given interac¬ 
tion to appear in the Hamiltonian it must not contain a 
rapidly oscillating factor ^ one finds that the Klein 

factors either cancel or combine to 

Taking this into account, we And five symmetry- 
allowed dimension two operators 

0^1 = COS$s , 

0 ^ 2 '’ = + {dxPs)'^ , 

of ^ = d^ifsdxPs , 

of^ = {d^^Pcf + {dxPcf , 

of ^ = d^ipc.dx'Pc ■ (A14) 

We note that none of these leads to a coupling between 
spin and charge sectors. In order to see that the operator 
cos$s is allowed, but sind)^ is not, one needs to consider 
the structure of Klein factors in (48). 


3. Dimension three operators allowed by symmetry 

The analogous analysis for dimension three operators 
gives the three possible perturbations involving only the 
charge sector 

of^ = {dxPcf + {dxPc)^ , 

of^ = {da:(pcfdoo(pc + {dxPc)'^ , 

of ^ = dxPcdlpc - dxPcdlpc, (A15) 

and four possible perturbations that couple spin and 
charge sectors together 

of^ = {dx^sfdxPc + {dxPsYdxPc , 
of^ = {dx^PsfdxPc + {dxPsfdx'Pc , 

Og ^ = dxPs^xPsi^xPc T ^xPc) : 

of ^ = {dxPc + dxPc) cos($s) . (A 16) 


4. Dimension four operators allowed by symmetry 

Finally we consider all symmetry-allowed dimension 
four operators. 


a. Charge sector only 

We And six perturbations involving only the charge 
sector 

of ^ = (dxPc)^ + {dxPc)^ , 

Of^ = {dxPcfidxPcf , 
of^ = {dxPcfdxPc + {dxPcTdxPc , 
of^ = {dlpcf + {dlpcf , 
of ^ = dlg^cdlpc , 

of^ = dlipc{.dxPcf - dlipc{dxPcf- (A17) 

Here we have used that 

dlPcdxPc + dl^pA^c = -2dltpcdlifc 

+ total derivative. (A18) 
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b. Spin sector only 

We find altogether eight symmetry allowed perturba¬ 
tions involving only the spin sector 

= {da,(ps)'^ + , 

oi^'> = , 

= {dx(Psfda;<fs + {dj;(fs)^da;(Ps , 

= (dlvsf + , 

o[f = , 

0 ^ 2 ^ = cos{'^s)dx^sdx(ps , 

0 ^3^ = cos($s) [{dx^sY + [dx^sf] ■ (A19) 

c. Terms coupling charge and spin 

Finally, there are eight symmetry-allowed perturba¬ 
tions involving both spin and charge sectors 

o[f = cos{^s)dxTcdxTc , 

= cos($s) [{dxTcf + {dxTcf] , 

+ {dxTsfWxTcf + {dxTc)\ 

= [{^xTsf + {dxTsY]dxTcdxTc, 

= dxTsdxTs[{dxTcf + {dxTcf'], 

^20 “ ^xTs^xTs ^xTc^xTc-, 

O2? = [{dxTsf + [dxTsf] [dlTc - dl(pc], 

0^22 = dxTsdxTs [d’i'Tc - (A20) 

Appendix B: List of irrelevant operators at the LE 
point 

The symmetries of the Hamiltonian at the Luther- 
Emery point are the same as those listed in Sec. A. 
Since at the LE point all symmetry-allowed operators 
in the Hamiltonian are local in terms of free holons and 
spinons, we shall give the list of irrelevant operators di¬ 
rectly in the fermionic representation. Among the op¬ 
erators listed in Appendix A, those which involve only 
derivatives of ip a, Ta preserve the same scaling dimension 
at the LE point; their expressions in the fermionic basis 
can be obtained straightforwardly by using bosonization 
identities for spinless fermions. 

On the other hand, operators that contain cos($s) re¬ 
quire a more careful analysis. First we note that, taking 
into account the Klein factors, the operator is actually 
represented by 

cos($s) cos($s). (Bl) 

Recall that spin flip and parity act nontrivially on 
However, the product of four Majorana fermions in Eq. 


(Bl) is invariant under both transformations. As a re¬ 
sult, the Klein factors can be safely omitted in the sym¬ 
metry analysis of the bosonized perturbations to the Lut- 
tinger model. At the LE point, however, cos($s) must be 
refermionized into free spinons according to Eq. (67). In 
terms of spinon operators, spin-flip symmetry is equiva¬ 
lent to a particle-hole symmetry 

Rg o RI, Lg o L\, (B2) 

while leaving the spinon Klein factors r]s,fjs invariant. 
Parity acts on spinon operators in the form 

PR,{x)P = Ls{-x), PLs{x)P = Rs{-x), (B3) 

Pr]sP = fjs, Pfj^P = r]s. (B4) 

The refermionization of cos(<i)s) at the LE point yields 

cos(4>s) —)• cos($*/-\/ 2 ) 

~ ’^{VsRIvsL, + r]sR,fj,Ll) 

= -Tr^TjsfjsRlL^+fjs'nsLlR^). (B5) 

Notice that the operator in Eq. (B5) is invariant un¬ 
der the spin-flip transformation (B2) only if we take the 
spinon Klein factors into account explicitly. Neverthe¬ 
less, in the following list of irrelevant operators we shall 
omit the Klein factors for short and write simply 

cos(<i>«)^(RlL, + LtRj. (B 6 ) 

Since all the operators that stem from cos($s) contain 
the same combination of Klein factors in Eq. (B5), we 
can adopt the prescription in Eq. (B 6 ) supplemented by 
the ad hoc rule that spin-flip symmetry takes R^ —>■ Rl 
but Lg —>• — 

Importantly, the operator cos(<i>s) has scaling dimen¬ 
sion 2 at the SU(2)-symmetric weak coupling regime, but 
dimension 1 at the LE point. Therefore, the scaling 
dimension of perturbations that contain cos(<i)s) is re¬ 
duced by 1 as we go from weak coupling to the LE point. 
For instance, the marginal operator 0\ ' = cos(<i)s) 
at the SU(2) point becomes the relevant mass term 
of the Thirring model, cos(<i)*/\/2) ~ R\Lg -|- L\Rg 
at the LE point, while the irrelevant (dimension-three) 

operator = {dxPc + dxPc) cos{^s) at the SU(2) 

point gives rise to the marginal spin-charge coupling 
(RlLg + LlRg){RlR^ + L\L^) at the LE point. This 
implies that, in order to have the complete list of irrele¬ 
vant operators up to dimension four at the LE point, we 
have to consider the refermionization of operators that 
have dimension 5 at the SU(2) point and were not in¬ 
cluded in the list in Appendix A. The latter correspond 
to operators through in the list below. 
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1. Dimension-three operators allowed by symmetry 


=RldlR, + LldlL^ + h.c., (B7) 

Vf ^ = R\R^L\id^L^ - L\L^Rlid^R^ + h.c.. (B8) 
= d,{R\RMW - BMlMRc, (B9) 
Vf ^ = d^Rld^L^ + d,L\d,R^, (BIO) 

Vf) = iR%R,LlR, - zL%lXLs + h-c., (Bll) 
Vf) = RlR.R^d.Rs - LlL^Llzd^L^ + h.c., (B12) 
Vf) = RlR^Llid.L^ - LlL^Rlid,R,+h.c., (B13) 
Vf) = {RlR, + L\lMRsL\L„ (B14) 

= (iR\dxRc - iLldxL^ + h.c.) 

xiRlL,+LlR,), (B15) 

= rIr^lIlMLs + LlRs), (B16) 

= d,{RlR, - lIlMl, + LlR,). (B17) 


2. Dimension-four operators allowed by symmetry 

a. Charge sector only 


vi") = iRldlR, - iLldlL, + h.c., (B18) 

= R%R^Lld,,L^ + h.c., (B19) 

Vf = RldlR.LlLc + LldlL^Rc + h.c., (B20) 

Vf ^ = R%RlR,d,R, + L%LlL^d,L^, (B21) 

vf (B22) 

Vf = d,{RlRM^d,Lc + dMLM^d^Rc 

-fh.c.. (B23) 


b. Spin sector only 


= iRldlR, - tLldlL,, (B24) 

Vf ^ = Rld^RMdxL, + h.c., (B25) 

= Rld^RlL.d^L^ + L%LlR,d,R,„ (B26) 

= RldlR^L, + lIoIlXRs + h.c., (B27) 

Vif = R%RlR,d,R, + Lld^LlL^d^L^, (B28) 

= d,{RlRjd,{LlL,). (B29) 



c. Terms coupling charge and spin 


y(4) 

*^13 

= {iR\dxRc + h.c.)(ii?Ji9a:i?s + h.c.) 



+(iL\dxL^ + h.c.)(iL\dxL^ + h.c.). 

(B30) 

y(4) 

'^14 

= {iR\dxRc + h.c.){iL\dxL g + h.c.) 



+{iL\dxL^ -T h.c.){iR\dxRs + h.c.), 

(B31) 

y(4) 

= R\R^L\L^(iR\dxRs ~ iL\dxL,, + h.c. 

),(B32) 

y(4) 

^16 

= (iRldxR^ - iLldxL^ + h.c.)RlR,LlL 

„(B33) 

V(4) 

|/l7 

= rIr^lIl^rIr^lIl^, 

(B34) 

y(4) 

= dxiRlR,)itRldxR,+h.c.) 



+dx{LlL^mldxL,+h.c.), 

(B35) 

y(4) 

^19 

= dxiRlRc)iiLldxL^ -f h.c.) 



+dx{LlLj{iRldxRs + h.c.). 

(B36) 

.^(4) 

^20 

= dx{RlR,-LlLjRlRsLlL,, 

(B37) 

y(4) 

'^21 

= {RldlR, + LldlL,+h.c.) 



^ (rIRs + rIRs)^ 

(B38) 

y(4) 

= {R\R^L\idxLc ~ R\RcR\tdxRc + h.c. 

) 


^ {R\Rs + R\Rs) 

(B39) 

y(4) 

= R^R^RlidxRsLlRg 



~R1RcRI'^^xLsRIl^ + h.c.. 

(B40) 

y(4) 
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= RlRcLUdxL,,RlL, 



-RlLcRltdxRsLlRs + h.c. 

(B41) 

y(4) 

= a|(i?ti?, + LtLj(i?tL^+Lti?J, 

(B42) 

y(4) 

= dxiRlidxR, + LlidxL^mlL^ + LlR^ 

,)(B43) 

^(4) 

1427 

= [dxiRlRMLc-dxiLlLMRc] 



^ (rIRs + rIRs)- 

(B44) 

Appendix C: Constraining the mobile impurity 


model in the SU(2)-symmetric case 

We can impose SU(2) symmetry in the parameters of 
the mobile impurity model of Section V by requiring that 
the edge exponents of longitudinal and transverse spin- 
spin correlations coincide'*^^. First consider the longitudi¬ 
nal component of the spin density operator: 

S^{x) ~ R\{x)R^{x) - rI{x)R^{x) 

r^0l{x)0,ix), (Cl) 

where the charge strings are cancelled in the sense of the 
lowest order in the OPE. We then project Os{x) so as to 
create a high-energy spinon: 

Os{x) ^ 


(C2) 
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while Ol{x) acts only in the low-energy subband: 


(C3) 

Cancelling the neutral string, we obtain 

S^{x) ~ (C4) 

In terms of the transformed impurity field 


(C5) 

Now, consider the transverse component: 

S-{x) ~ Rl{x)R^ix) 

~ Osix)Osix) 

_ (C 6 ) 


Given the expressions (C5) and (C 6 ) we can calcu¬ 
late threshold exponents in the Fourier transforms of 
the spin correlations functions (^^(x, t)S'^(0,0)) and 
(5'+(a;,t)S'“(0, 0)) and impose that they share the same 
exponents in the SU(2)-symmetric case. This has to be 
the case even for higher harmonics, taking into account 
backscattering processes^^. A shortcut is to compare the 
scaling dimensions of the strings in Eqs. (C5) and (C 6 ). 
We must have 


It follows that SU(2) symmetry imposes 


(C7) 

(C 8 ) 


Is = Is = - 


1 

4^2' 


(C9) 


These values of 7 s, 7 s are of order 1, as expected since 
the spinons are strongly interacting at the SU(2) point. 
Moreover, these values are consistent with the exact re¬ 
sult for the Hubbard model at zero magnetic field (see 
Sec. VB). Using Eq. (C9) and setting = 1 in Eq. 
( 110 ) yields 


z/s,+ = Vs,- = 0. (CIO) 


Appendix D: Priedel oscillations 

In this appendix we summarize some useful facts re¬ 
garding Friedel oscillations for open boundary conditions. 
At low energies the charge density operator has an ex¬ 
pansion of the form 

•nj ^^P2nkF{x), (Dl) 

n =0 


where p 2 nkF (a:) denotes the Fourier components with mo¬ 
menta ±2nkp. For periodic boundary conditions it fol¬ 
lows from (48) and the site-parity symmetry (A9) that 
the leading contributions to the 2 kF and Akp components 
are 


P 2 kF{x) = A 2 kF sin 





Pikp {x) = cos (4>c 


^kpx) -I- ..., 


(D2) 


where A 2 kp and ^ 4 ^^ are non-universal amplitudes that 
include Klein factors. The leading contributions to zero 
momentum component are 


Po{x) = n- + Aod^<^c 

ZTT 



(D3) 


where n is the band filling. 


1. Open boundary conditions 


For open boundary conditions we still have expansions 
like (48), but the chiral Bose fields no longer commute 
and hence the Klein factors F^ ^ need to be adjusted 
accordingly. The mode expansions for the chiral Bose 
fields are 


Vc{x) 


(Pc{x) 


a + TTo 
2 
00 

n—l 

a - TTo 
2 
00 

-E 


X(po 

2{L + 1) 

I - h.c.] , 

xipo 

2{L + 1) 

I - h.c.] , 


(D4) 


where [ttq, 1 ^ 0 ] = Sttz and = pnjiL 1). The commu¬ 
tator between different chiralities is 


{ 0 if X = ?/ = 0, 

2pi if 0 < X, 2/< L-I-1, (D5) 

47 rf \i X = y = L + 1. 


The boundary conditions on the Fermi creation and an¬ 
nihilation operators are 


C 4 — 0 ,cr 9 (^9) 

Imposing the boundary conditions (D 6 ) on the bosonized 
expression (48) (with Klein factors appropriate for the 
mode expansions (D4)) leads to conditions of the form 

®c(0)|V^o) = CilV'o) , ^c(7- + l)l^o) = C2|'0o) , 

$,(0)|V'o) = ^s{L + l)|V'o) = cMo) ■ (D7) 

Importantly, the actual values of the ci _2 and c] 2 depend 
on all the amplitudes A„ in (48) and are thus non- 
universal. 









26 


2. Friedel Oscillations 

For open boundary conditions the 2k f and Akp com¬ 
ponents of the total charge density are again given by 
expressions of the form (D2). Using the mode expansions 
we then can determine the form of the Friedel oscillations 

sin (2k'pX — 

{'ipo\p2k,,{x)\'ipo) ~ B2kp ’ 
{^o\p4kir{xMo) ~ ^2K ’ (D8) 


where e is a short-distance cutoff, we have taken into 
account Kc ^ 1 and defined 


kp — kp -\- 


Cl - C2 
4{L + 1) ■ 


(D9) 


The upshot is that Friedel oscillations involve two non- 
universal interaction-dependent parameters: the phase 
shift Cl and the 1/L shift oikp. 


^ C. Xu and S. Sachdev, Phys. Rev. Lett. 105, 057201 

( 2010 ). 

^ V. V. Deshpande, M. Bockrath, L. I. Glazman, and A. 

Yacoby, Nature 464, 209 (2010). 

® O. M. Auslaender et al., Science 295, 825 (2002). 

^ B. J. Kim et al, Nat. Phys. 2, 397 (2006). 

® Y. Jompol et al., Science 31, 597 (2009). 

® J. Schlappa et al., Nature 485, 82 (2012). 

^ F. H. L. Essler, H. Frahm, F. Gohmann, A. Kliimper, 
and V. E. Korepin, The One Dimensional Hubbard Model 
(Cambridge University Press, Cambridge, 2005). 

® F. H. L. Essler and V. E. Korepin, Phys. Rev. Lett. 72, 
908 (1994). 

® F. H. L. Essler and V. E. Korepin, Nucl. Phys. B426, 505 
(1994). 

N. Andrei, arXiv:cond-mat/9408101. 

T. Deguchi, F. H. L. Essler, F. Gohmann, V.E. Korepin, A. 
Kliimper, and K. Kusakabe, Phys. Rep. 331, 197 (2000). 

K. Penc, F. Mila, and H. Shiba, Phys. Rev. Lett. 75, 894- 
897 (1995). 

J. Favand, S. Haas, K. Penc, F. Mila, and E. Dagotto, 
Phys. Rev. B 55, R4859 (1997); 

K. Penc, K. Hallberg, F. Mila, and H. Shiba, Phys. Rev. 
B 55, 15475 (1997). 

H. Benthien, F. Gebhard, and E. Jeckelmann, Phys. Rev. 
Lett. 92, 256401 (2004). 

A.E. Feiguin and D.A. Huse, Phys. Rev. B79, 100507(R) 
(2009). 

M. Arikawa, Y. Saiga, and Y. Kuramoto, Phys. Rev. Lett. 
86, 3096 (2001). 

M. Arikawa, T. Yamamoto, Y. Saiga, and Y. Kuramoto, 
Nucl. Phys. B702, 380 (2004). 

A. B. Zamolodchikov and Al. B. Zamolodchikov, Ann. 
Phys. 120, 253 (1979). 

L. D. Faddeev, Sov. Sci. Rev. Math. Phys. C 1, 107 (1980). 
R. A. J. van Elburg and K. Schoutens, Phys. Rev. B58, 
15704 (1998). 

G. 1. Japaridze, A. A. Nersesyan, and P. B. Wiegmann, 
Nucl. Phys. B230, 511 (1984). 

F. Woynarovich, J. Phys. A 22, 4243 (1989). 

H. Frahm and V. E. Korepin, Phys. Rev. B 42, 10553 
(1990). 

I. Affleck, in Fields, Strings and Critical Phenomena, eds 
E. Brezin and J. Zinn-Justin, (Elsevier, Amsterdam, 1989). 
A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, 
Bosonization and Strongly Correlated Systems (Cambridge 


University Press, Cambridge, 1999). 

T. Giamarchi, Quantum Physics in One Dimension (Glare- 
don Press, Oxford, 2004). 

M. Pustilnik, M. Khodas, A. Kamenev, and L.I. Glazman, 
Phys. Rev. Lett. 96, 196405 (2006). 

A. Rozhkov, Eur. Phys. J. B 47, 193 (2005). 

A. V. Rozhkov, Phys. Rev. B 74, 245123 (2006). 

J. M. P. Carmelo, K. Penc, L.M. Martelo, P. D. Sacra¬ 
mento, J. M. B. Lopes Dos Santos, R. Claessen, M. Sing 
and U. Schwingenschlogl, Europhys. Lett. 67, 233 (2004). 
J. M. P. Carmelo, K. Penc, P. D. Sacramento, M. Sing and 
R. Claessen, J. Phys.: Condens. Matter 18, 5191 (2006). 
J. M. P. Carmelo, D. Bozi and K. Penc, J. Phys. Cond. 
Mat. 20, 415103 (2008). 

R. G. Pereira, J. Sirker, J.-S. Caux, R. Hagemans, J. M. 
Maillet, S. R. White, and 1. Affleck, Phys. Rev. Lett. 96, 
257202 (2006). 

E. Bettelheim, A. G. Abanov, and P. Wiegmann, Phys. 
Rev. Lett. 97, 246401 (2006). 

M. Khodas, M. Pustilnik, A. Kamenev, and L.I. Glazman, 
Phys. Rev. Lett. 99, 110405 (2007). 

M. Khodas, M. Pustilnik, A. Kamenev, and L.I. Glazman, 
Phys. Rev. B 76, 155402 (2007). 

R. G. Pereira, S. R. White, and 1. Affleck, Phys. Rev. Lett. 
100, 027206 (2008). 

A. Imambekov and L.I. Glazman, Phys. Rev. Lett. 100, 
206805 (2008). 

V. V. Cheianov and M. Pustilnik, Phys. Rev. Lett. 100, 
126403 (2008). 

E. Bettelheim, A. G. Abanov, and P. Wiegmann, J. Phys. 
A41, 392003 (2008). 

A. G. Abanov, E. Bettelheim, and P. Wiegmann, J. Phys. 
A42, 135201 (2009). 

A. Imambekov and L. I. Glazman, Phys.Rev. Lett. 102, 
126405 (2009). 

A. Imambekov and L. I. Glazman, Science 323, 228 (2009). 
R. G. Pereira, S. R. White, and 1. Affleck, Phys. Rev. B 79, 
165113 (2009). 

R. G. Pereira, Int. J. Mod. Phys. B 26, 1244008 (2012). 
A. V. Rozhkov, Phys. Rev. Lett. 112, 106403 (2014). 

T. Price and A. Lamacraft, Phys. Rev. B90, 241415 (2014). 
R. G. Pereira and E. Sela, Phys. Rev. B 82, 115324 (2010). 
R. G. Pereira, K. Penc, S. R. White, P. D. Sacramento, 
J. M. P. Carmelo, Phys. Rev. B 85, 165132 (2012). 

T. L. Schmidt, A. Imambekov, and L. 1. Glazman, Phys. 
Rev. Lett. 104, 116403 (2010). 






27 


T. L. Schmidt, A. Imambekov, and L. I. Glazman, Phys. 
Rev. B 82 245104 (2010). 

F. H. L. Essler, Phys. Rev. B 81, 205120 (2010). 

A. Imambekov, T. L. Schmidt, and L. I. Glazman, Rev. 
Mod. Phys 84, 1253 (2012). 

O. Tsyplyatyev and A. J. Schofield, Phys. Rev. B 90, 
014309 (2014). 

L. Seabra, F. H. L. Essler, F. Pollmann, I. Schneider, and 
T. Veness, Phys. Rev. B 90, 245127 (2014). 

V. E. Korepin, N. M. Bogoliubov, and A. G. Izergin, Quan¬ 
tum Inverse Scattering Method and Correlation Functions 
(Gambridge University Press, Gambridge, 1999). 

L. D. Faddeev and L. Takhtajan, Jour. Sov. Math. 24, 241 
(1984). 

F. H. L. Essler and R. M. Konik, in From Fields to Strings: 
Cireumnavigating Theoretical Physics, ed. M. Shifman, A. 
Vainshtein and J. Wheater (World Scientific, Singapore, 
2005) and references therein; cond-mat/0412421. 

K. D. Schotte and U. Schotte, Phys. Rev. 182, 479 (1969). 
A. Luther and V. J. Emery, Phys. Rev. Lett. 33, 589 
(1974). 

F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1982). 

H. Frahm and V. E. Korepin, Phys. Rev. B 43, 5653 (1991). 
D. B. Abraham, F. H. L. Essler, and A. Maciolek, Phys. 
Rev. Lett. 98, 170602 (2007). 


S. Coleman, Phys. Rev. Dll, 2088, (1975). 

M. B. Zvonarev, V. V. Cheianov, and T. Giamarchi, J. 
Stat. Mech. P07035 (2009). 

H. Karimi and 1. Affleck, Phys. Rev. B84, 174420 (2011). 
C.G. Callan, I.R. Klebanov, A.W.W. Ludwig and J.M. 
Maldacena, Nucl.Phys. B422, 417 (1994). 

T. Xiang, Phys. Rev. B 53, R10445(R) (1996). 

U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). 

F. H. L. Essler, V. E. Korepin, and K. Schoutens, Phys. 
Rev. Lett. 67, 3848 (1991). 

F. H. L. Essler, V. E. Korepin, and K. Schoutens, Nucl. 
Phys. B372, 559 (1992). 

M. Ogata and H. Shiba, Phys. Rev. B 41, 2326 (1990). 

A. W. W. Ludwig and K. J. Wiese, Nucl. Phys. B 661, 
577 (2003). 

S. Eggert, Phys. Rev. B 54, R9612 (1996). 

K. Okamoto and K. Nomura, Phys. Lett. A169, 433 
(1992). 

S.A. Soiling, M. Bortz, 1. Schneider, A. Struck, M. Fleis- 
chhauer, and S. Eggert, Phys. Rev. B 79, 195114 (2009). 
S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

S. R. White, Phys. Rev. B 48, 10345 (1993). 



